{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "sirc! (generic function with 1 method)"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "include(\"SIRcasedeath.jl\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " 37.554100 seconds (262.95 M allocations: 12.673 GiB, 8.53% gc time)\n"
     ]
    }
   ],
   "source": [
    "csvfile = CSV.File(\"covid_timeseries_UnitedKingdom.csv\")\n",
    "df = DataFrame(csvfile)\n",
    "df[:,:active] = df[:,:cum_cases] - df[:,:cum_recover] - df[:,:cum_deaths]\n",
    "# @time out = mcmc(dfuk,[8,3,4],ones(6),6.67e7,10);\n",
    "@time out = mcmc(df,[8,3,4],ones(6),6.67e7,100000);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAGdCAYAAADJ6dNTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdd3gUVRcG8Hc3PaTQQ0novUmH0KRIESkiCEgXpEmR9qk0BRQDKlUwIC2C9N4EKdJ7C70jBEhCDVkS0ne/PybbsrMt2Zq8v+fJkyl3Zk5E2JM7954rUSgUChARERE5CKm9AyAiIiLSxOSEiIiIHAqTEyIiInIoTE6IiIjIoTA5ISIiIofC5ISIiIgcCpMTIiIicihMToiIiMihuNo7AFPI5XJERkbC19cXEonE3uEQERGRCRQKBd6+fYsiRYpAKjW9P8QpkpPIyEgEBQXZOwwiIiLKhMePHyMwMNDk9k6RnPj6+gIQfjg/Pz87R0NERESmkMlkCAoKUn2Om8opkhPlqxw/Pz8mJ0RERE7G3CEZHBBLREREDoXJCRERETkUJidERETkUJicEBERkUNhckJEREQOhckJERERORQmJ0RERORQmJwQERGRQ2FyQkRERA6FyQkRERE5FCYnRERE5FCYnBAREZFDYXJCREQ5nkKhwOX9f+PJrev2DoXgJKsSExERWdOjq+E4sPR3AMDY9bvsHI1tJL2Lh1wuh5ePr71D0cGeEyIiyhHevn6Jh+EXoFAooJDLAQCxz6Nx/8IZ7JozQ9Xu+cMH9grRZhRyORZ83g2/D/gMKclJ9g5HB3tOiIgoR/hjaD+t/S9+W4qlI77Qabfqm5H4fM5i5C1SFMmJCTi5cQ3K1WuIIuUq2ChS60tLTVVtx79+jdyFCtsxGl3sOSEiIptLSU7CjWOH8E4Wa9Z1MVFPsf+PBXjzLDrLMaz97mu9556mjz3ZNWcGLuzairWTx4m2u7hnB2Z1a4cT61fpnHsni4Xs5Yssx2ltca9f2TsEHUxOiIjI5o7+tQJ7FszCxh8mmnXd+inf4srBvdgS8l2WY4iPea33nKuHBwDgv/ALetvIXr7AobA/AACnt6xXHb975iT2LJiF0IE9sWTY54iLeY3U5GSs/GYk/l2xOMtxZ4ZCocCZrRuwdeZUpCQlQgGF6ty2X3/AjWOHEPv8mV1iE8PkhIiIbO7O6eMAgJcRD826Lv5NDAAgJioSKYmJlg5LxdXdHQqFwmAb2cvnWvunNq+FQqHAjtk/4caxQ6rji4f0wbzen+DFwwe4tHenVeI15uTGNTi+biUeXDyHi3t2IikuTnUuKT4eexbMQuTdW3aJTQyTEyIiyrTbp45hzeRxSHgrM6l9TNRTXD20D/L0Aali0lJTsXriGOz74zfVsVdPHuPSP9qzaHbN/9mkZ0Zcu4xTm9aa1FZpx6/TcfHv7QbbXNqjnWic3LAas7u3N3pvZQ9FTHQklo8ahLPbN2mdT4yLw4LPu2HRkD5ITkwwK259Tm9W//zH1/6JiGuXddq4uDrOMFQmJ0RElGm75s5E1J1b+GfRfNHzD8MvYP+SBUhJEno5lo8ajH2L5iPRQDJz79wpRN+7g6sH/1EdCxs7FP8uX6TV7sGFs8L3i+fw57hheHLjmuj9Nv4wESc3rjbr5wKAwyuXau1HXLuCB5fOAQCePbin6v0xl7L3Z9+i+YiJisSxNWE4sHShqqfm8MqlSHoXj/iY1zoJkKW4uLnrHJO6OE5y4jiREBGR07p//jTObt+Euh27qI6d3b4Jx9aEAQBc3NzQvN9gk+514+i/Zj173+L5iH8Tg31//Ib+c9VjOvYvWYArB/aadS9DNv4wAQBQonotPDQwFsWYlMREJCcm4MlNdTJ1ef8elKkTDDcPT1w/ckB1PDnhHQAgNSUFLq6ukEgkZj/vhcirM29/f51jLi4uZt/bWpicEBGRRRxbEwa//AXgky8//AsEqBITQHgF8vrpE9HrUpIS4ebhqdqXvVCP5VDI5ZBIDXfyK3si3maYGWPJxERTVhITAIi4Fo5N0yfpHN/8k+4g3+gH95D0Lh7LvhqEwmXKodM335v9vJX/G65zTCHyWk3539ERMDkhIrIjhUKB5w8fIG+Rolof0NaWEPcW6777GuWDG6HBpz0tdt/d83/Re+7RlUuix28cPQQ3T0+Uq98Irm5uePn4keqcXJ4GeXKyWTGkJCbi1ZMIs66xJc0aI8ZIJBKEDuyJtNRUPLh4zqoxGEsCbYnJCRGRHd06cQR///YrCpUph57TZ9vsuWe3bcTrp49xatNa1OvUFe9iY+GbL79Z93j24J5FYjmwdCEA4NWTCDT+rK/Wub/n/4o7Z07ovfbIX8tV26kpyXhw8RwOLl8E2QvHmRab0YXd20xumz+omN6kTszryCdYMXoIAKBWu06o/dHHou2e3rqhcyxP4aImP8faJApjc6UcgEwmg7+/P2JjY+Hn52fvcIiILGb5qEGIiYoEIL6mS0zUUyS9e4dCpcta9LlzenwMeZrw23OB4iXx4tF/6DVjHgJKljbperk8DSv/N8KiPRR5igSiwac9sHueabNwcqKx63fh6a0b8C8YAJ+8+XTOz+rWLtP3HrNuZ6bGtBiS2c9vx+nDISLKgdy9vA2eXz5qMFZPGI1zO7dgVrd2iL5/1yLPVSYmAPDi0X8AgPB/TFvw7tmDe1jY/zOLvzqJiXzCxMSIJzevYd33X2Px0L549eSxRe9t6cQkK5icEBFZ0evIp9g5OwTP/rsvet5QvQ9NR9NfX6yeMFp1LP5NDC7u2YnTW9Yj6u5tnWvEBj0a4uruoffcswf3cGrTWqQmJ2Ptd/9TzSLJTmp+2AHDlq+zdxgGrZ/yrWo7bOxQnfM+Zr6aU2rUvU+mY7IGjjkhIsoihUKB108fI0/hopCmT8e8feo4DoUtVs2AuHPmBHr8OAsBpcrg5eNH8Pbzh0/efHihsQLuma0bULhsBcjlabiways6jJkg+rzE+Dh45vLBspEDVfVDTqxfhbHrdyE1ORnH161UjWuo37k7GnbtpXV90jvxxOKdgdkaf40fBQCZqhfiaPr+uhB/jhumdeyT8VNRvGp11Z+fM3l8/QpePX2C9z5og7hXLzN1j8rvt7BwVFljVnISGhqK0NBQPHz4EABQuXJlfPfdd/jwww9F24eFheHzzz/XOZ6QkABPT9uNSicisqYjq5bhwu5tqNSkOT4cNgYAsGvuDJ12ayaNReEy5RF1T+jlGBG2Qev88XUrtfbn9+0CMQv7d0elxs1UiYnSmW0bcXztn1rHTm9ep5WcpCQlYsHnXUXve+fMCSjkckRcvwIXV1cEVqyC15FP8ejKRdH2zsrTxxdj1u3E5f17cHDZ7wCAktVr2TmqzJndvQMUCqGH7O6ZzBWFq9aijej4FXsy67VOYGAgZsyYgfPnz+P8+fNo3rw5OnbsiOvXr+u9xs/PD1FRUVpfTEyIKDu5+u8+AELxsPsXzhp8naJMTADgefpYj8zQXLtFKWNiIiY+xnAti/n9PsWmHydh/ZRvkZzwDitGD7bbYnXWIpVKIZFIkDugkOj5Wh91NPueDbv1zmpYKvU6dTO5rTIxAYQKtpnRcpBuHRR7Mys5ad++Pdq2bYty5cqhXLlymD59Onx8fHD69Gm910gkEhQqVEjri4goO9FMRrb9PA1XD+036bqbR3UTDGtLjI8zeD41KUm1/Vs/8R4Wp5c+8LN4tRpo3n8I2o36Vuv0+72/gIubm1m3rP9JN1Ro+L5FwmvU3XKJjrPK9IDYtLQ0rFu3DvHx8QgODtbbLi4uDsWLF0dgYCDatWuHS5eMz9dOSkqCTCbT+iIiclQZe0r2ayxYZ8iVg9apYJqRPC0N9y+cxcU9O3Bq0xqbPNNaStaojS6TfrTIvSQSCWq0bofywY10jn808n8m36di42YAAA9v3ZlXysX0PH390LhHv8wHm0X5i5Ww27Mzw+wBsVevXkVwcDASExPh4+ODrVu3olKlSqJtK1SogLCwMFStWhUymQzz5s1Dw4YNcfnyZZQtq3/OfkhICKZOnWpuaEREdiGXp9k7BINePPoP236eZpdnf/bDr4iPjcGOX6db5H4dx02yyeq5ZWrXR4027VG4bHnEPn+GE+tXoc8vC7Br7ky8fqo9hVc5zsg3XwHVscpNP0Dk7ZvoFTIHVw7+g3L1GuLRtXCjzy1hpbEvfWbOR/SDu1gzcaxV7m9pZv8Jly9fHuHh4Xjz5g02b96Mvn374siRI6IJSv369VG/fn3VfsOGDVGzZk389ttvmD9ffAVLABg/fjzGjBmj2pfJZAgKCjI3VCIiq3n5+BEO/bkEDbr0QNHylfD4xlV7h6TXhb+32+3ZCrkcJavXtsi9qjZvZZPEBBBKuTf/XL1QYf1PhHEgfX7+DZf37cahP5eo26a/JnqvVVu8evoYpWrWQYUGTVTna7frBACQStUzgWq164QLu7bqPLdZ30E6x4pVqZbp8SSaP0/hMuWzdA9bMvu1jru7O8qUKYPatWsjJCQE7733HubNm2faw6RS1KlTB3fvGi4i5OHhAT8/P60vIqKsePMsGhumTcDDy5aZebLt52mIuBqOdd9/7dCJCQDcPa2//Lu1KRRyuGYYv6FcuTi4Sw+zBpI2769b18Mc/gUD4F8wAF4+vpm+h4urK6p+0AYeuXLpnPPM5YO2w8dqJSaaygc3RmClKqjYqCkadhVfz8g3fdaMb351L8ynk3/KdLwZVW+d+QqytpTlFFShUCBJYwCVsbbh4eGoWrVqVh9LRGSy1ORkLBv5BQChJoRYmXhTKRQKSCQSyDKsgOvIUlPMWzjPkvIFFgMgjHl4GfEQANDos76o2qIN/AsGIC01FSfWrzLpXppJTv6g4loLBJqi/9w/AGR9gTs3dw8MX74ej29cRa7ceU2+ztXdHd2+151irnXv9NmsDbv2wt7f52QpTqXAilVU242694FCLsfl/X9b5N7WYtaf0IQJE3Ds2DE8fPgQV69excSJE3H48GH07ClkgH369MH48eNV7adOnYp//vkHDx48QHh4OAYMGIDw8HAMGTLEsj8FEZEBp7es19q/e/YkLuzeBuXSYgq5HAqFAk9uXENKYqLYLQAIK7mu+noEdsz+CU6wLJlDUM568fBW9zQop/FKJBK4urmh/7w/RK/1yaP/g9/VQ38124zyFC6K9mPGQ+riYtEia0GVqiJvEessluflq/3GYMC8JWg1ZGSm7tV+jPpz2cPbGx988WWWYrMFs5KTZ8+eoXfv3ihfvjxatGiBM2fOYO/evWjZsiUAICIiAlFRUar2b968waBBg1CxYkW0atUKT58+xdGjR1G3bl3L/hRERHqkJCbizFbt5GTHrJ9weOVS3Dp+GE9uXcfszzpgdvf2WD/1W8zv20XvmjFRd27hRcRD3D1zErBhcjJ8xXq0HvKV6Lmu3/2EsnUb2CwWQD0AVJPmB6AmVzd3AEDLgcPg7Z8bzfoN1mmTp1ARjFm3E31+WYC6H38KAKj2QRsMXrRSp61S0fLqcY6BlaqgUJlyqv2hS9RVbN9r9RH6z12McvUaGvmpbK/loBEoU6c+Bi9aiSLlKqLNl+qlCUpWr4UKDd9Hk55CIdPchQqjarNWmXqOt5+/zjHln4uj4qrERJStndy4Ru/0WXcvL+QtEii6mN6YtTuQnJiAh5cvolTNOnDz8MTR1Stwbsdma4esQ/kaSmzF2ZYDh6NKs5Y4uno5Luw2b+BrxcbNcDO9mFvxajXw6IrxUg/KeDRj6TZlBgIrVsGNo/9iz8LZorED6ldihqQmJyPi+mUEVqgMdy9vzO3VCWkpKTr3SklMxPy+XeDq5o4vl67Bk1vXsSXkexSvVgNdJv6AB5fO4dqh/fjgi2GiH87OytRVh5t/PhhPb99E3Y5dULBEKZ3zUXdv4+DyRWjaZ4DWax9Ly+znN9fWIaJsLS7mld5zyQkJelf5nf1ZB9V25fdboEHXXplKTHIHFMabZ1HGG2ZSQKkykLq4oHqrdlrJSdXmrVSVa/XRTBS6TPzBpA++z+csAgAUKV8JkbdvoMPYCaoPt0pNmmslJ/4FA/Q+Tx9Xd3eUqlFHta9MTDJy8/TUSlZKVq+FQb+HwdtfSERK1aijdZ/sou+vCxH7/BlSEhOwe/4vom1cPTxQo0171GjTXu99Cpctj14hlhnTYg1clZiIyIjrRw7ixtF/M3VthYbiMzcM+WS86XWe3L28AAjd/p2+/V513MWEbns3M8ZtKOUtEggA6DBmPDpPmIYytetrna/avBVyBxTGiD834ovflpl9/4yKVa1uclvffPnh4mpeZVdnkz+oOErXqqu1mnX/udrLC2Rc6NEZseeEiBzerRNH8OLRfyhQvCTyBRbDkb+Wo1DpcpC6SBHcpYfB38gVcsu8uTZlRknTPgNxeOUSrWM1PuygNSC3Rf+hOLg8VKtN2boNcPfsSfjkzYcB85aYVTo9T2H1gEzNngJDg0kBYPTa7UiQyRB5+yaqtmht8vOUcuXOg1y58+gcbzU4c4M29Wk7fCzObtuIah+ILzCbU2lWJc5TuCjyBRZTjZVyxpWVM2JyQkQOQ6FQIOruLeQLLK4qBf5OFivafa0cHyEkLMVFZ03I09Jw7ZDhVxuWVKZOfZ3kJON4h+qtP4KHjw/+Tv+Z+vyyAAWKlUBC3Ft45vLRSbSqt/5ItV2h4fu4deKIar+AyFiCj776GvfPn0HNjzrqrHKs1HbEOEilLsiVOw/6/LJAdTygVBk8e3DPxJ/WNnLlzoNm/XQLk+V0hUprV1nv++tCzO4uvMbxzEIdF0fB5ISILCYlKRGvI5+iYIlSRscXJL2Lx4LPu6Fio6ZoO2IcAKGH5O/ffgUAjAjbAHcvb6z7zvAaJztmCQWq2o8Zj3L1GiIxPg7uXl6Qp8lxYOlCC/xUpvMvGIDGPfrh2Jowg+0qNnwf3n7+cHVzR4H0NU8yFgYbs3YHUlNT4OaufvXS5svRKPFeTVX9i0qNmurcu0KDJnqLgCmlJIlPl24xYKhOefN6nbrizNYNBu9HtpcvsBh6Tp+NXOk9ZBKJBC36D8WTW9eN/vk7AyYnRGQxm36cjMg7N9Fu1Lc6i6klvJXh9qnjuPj3drQb9Q2Orl4BALh5/DBaD/0KLq5uOLD0d1X73/p1RdM+XyAmKtKkZ5/bvgkBJUtj6Ygv4JnLx+jqu9ZSsXFTneSk/ehvsXPODLTQqHBa3MhYColUqpWYAEJ10srvt0DeIoF4cOk8qhsY8GiIvkGmBYpr98RUadYKwV16qJITzV4csj/N6dOA8OeTXf6MmJwQkcVE3rkJADi47Het5CQtNRW/f9FDtb973s9I1fiAPPrXCjTrNwjJCe+07nd45VKTnx19/y7O7dgCAHZLTACIjsMoV78Rvlq1Ba7ulqktUbhseRQua/46KR65ciEpPl7v1FGXDGMVWqcX/Rq1eiui7t1BYIXK5gdLlAmcrUNEFpfwVoa9oXNV1VYzVl19HfkEshfPVPsX9+zAyY3itUjMYawkt0/6uiXWJJW6oN2obwAAHcZMUB23VGKSFZ9N+wWf/fArChQvKXpes6x7A421X1xc3ZiYkE0xOSEiq7h++ADO7RTqgsjTUo2211cozZJ6z5iHgFJljTfMovLBjTFm7Q6UrWfbyq0ZZSxTni+wGIqUq2DStfU7dbNGSEQm4WsdIsq04+tWIjHuLZr1G4S0VN0E5O0roQCaPReeU/LJmw/e/rlRrEo1PHtgeGV0S8jq4nKW4OljfkXtrCyKSGQpTE6IKFMOLl+E8H+ED7L8QSVwbG2YThvlhJ2URNNWLrcm5TgLfVVee/w4CxumTUBqsv1jtRx1jZdcRuqeEDkS+6f2ROSQou7dxrrvv0HUvdui55WJCQDcO38ayQkJOm2u/rsPDy6dg0KeZrU4TRXcpYfB84XLlsfg0D/RqHufTN2//ehvM3WdNWlWif1ivumDi4nsjT0nRCRqzaRxgEKBNRPHYsjiVaKzUJQSZDK957bOmIpeM+ZZI0SzmLK0vaePD+p16ooi5Sog8V08dvw6XbSdu5cXUpKS0OOHX+GRKxeiH9xDWQdc9bbEezVRrEo1FCxZxiEG5BKZiqsSE5GojIvAtRoyErn88yDiWjia9OyPOT06mnwviUQKhUJuvKEVKcdSXD9yUFXETKlh116o37m7zjUKhUJVdRMAOv5vMhLexqJcvUaqCrZEpB9XJSaiLElLTUXknZsoXKa86G/Z+xbNV23nLmS8F0KTvRMTTWXrBuskJyVr1BZtq1nltk6HzihTu55VYyMiAcecEBlw99wpnNu5xd5h2MSxNSuwYep4/LPI+CuYKwf22CAibV0m/oiydYWpuVIXV3SbOhMD5i9Fy0HDzbqPu5c32o4YhxofqntETFkorVTNOkbbEJFlsOeEyADlmIPACpXNrsj59vVLHF+7Eg279oJfgYLWCM+iLuzeDkBY3+aDL4YZbPvi0X+2CElL8WrVEVipMk5tWocS1WuqioK9eKi/q9jFzQ3efrl1BqtWbNQUbp5euLRnJwDTkpPchYpkIXoiMgeTEyI94t/EqLZjn0ebnZz8MbQfAODG0X+drnbEgs+7Wu3eHf83Gdt/+SFT17q4uqFR995ax1JTddeJqd+5OxLj4tCi/xC993oXq/7zlafpn03Ub9bvSHr3Dj6ciktkM0xOiPTQHCsuVmDMkLiY15YOJ9swddxGgWIl8CLiIQCgwac99bZLfKs9UyhfYDE07NrL6P3dvdQDWj1y5dLbLl9gMaP3IiLLYnJCDis1JQWxz6ORr2iQXZ6vWUVUX3IS+zwapzatQ+32nZA/qDheRDzEma0bEJ2hNsjDyxdx5/RxNPqsL7z9/K0atyP7cNgYnWMSqRQKuTBgduz6XapZQr4FCuLjb75H7PNoBFWqqveemklknQ6dTS4ZL9Wo4OqX3/FfuxHlJExOyGFtnTEFEdcuo8OYCTZdo+SdLBZ7FszCw8sXVccKFCsh2nbrzGl49SQCd86cwMg/N2Ll/8QHZ27+6TsAQGJcHDqMnSDaJjvwyZsPg0P/xMU9O3EobLHOebFXY3kKFcHryCeq/bYjxuHC7u1o3m8w/PIXgF/+AgafqVkmvknPz02OVSLhfAAiR8W/neSwIq5dBgDsmP2TWddF3buNnbNDEPs8Wuv4gWWhWD/1W6QkJeLt65ei1ya9i0fowJ5aiQkgrJqrSaFQ4Nqh/Xj1JAIAkJKYYFKcyvbZVe+ZwnRjfQNM8xTWnYJcvfVHAAC/AgEAhMGqvULmwL9ggEnPrNCgCVzd3FHazGm+JWrUQu6AwigX3Nis64jI+thzQpmm7E7XrAWhz7KvBuJNdBRGr90OqVT44IqJegqffPnh5u6R6RiSE97h/K5tKB/cCPkCiyEtNQVrJo4FIMyW6fHjLKQkJmJ+3y6qa+b3Ebb7/PybztLx+lbGvXXiCD4a+T/V/qLBvfEu9o1Wm7tnThqNV7PioUKhwMmNqxH/JgYtBw43+t8xMT4Onrl8jD5D9uI5Di4PRe12nRBUuZrR9gBwZusGk9oZo3xldePoQZPat+g/FNVbf4Q8RQJRsESpTD3Ty9cPw8PWQ+pi3j9nbu4e6D93sXoBICJyGOw5oUxRKBRYPWE0/hj2uWq8gD7ytDS8iY4CAJzdtgmA8EG+fNRgzO/dGYfC/sCbZ9GGboEDSxfideQTzOrWDnN6fKw6fnT1CpzatAZhY7+EPC0NC/qpl3lXPnP91PGi97z67z6dYzFRkQZ/DqWMiYmp4jUGyu6aMwOnN6/D1YP/4MmNqwavO7HhLyzs3x13zpww+owlw/vjwcVz6YvYGV8N+NmDezi+bqXBNm2+HG30Ppqi7oqvx6NUo0175A4ojIqNmwIASlSrkaWxOC6ubiYlyRlJpNJMXUdE1sXkhDIlNSUZzx7cQ9yrl7h54ojedpF3bmktEHdi/SqkJidrvQK5uGcHlo38wmCSc3n/HqwYLUwLlaepB6dG3r6p2t7440Skpqg/jBPSZ3FoDmzV5ObpqXPswcVzemNQvmbSHB9hruSEd6ptzUTj1ZPHBq87vXkdAODgslCznjev9yd4cvMaLu/fg39XLIZCLkf4P7sxq1s7PLl5DQDw1/hRRu9T+f0WZj3Xy0ii0fzzweg/7w94eOufJUNEOReTE8qUFw/VRbgu7d0p2ib2+TOsnTwOh/5conX8zLaNou2j74snEWIeXQkX4kifagoAT25cM/l6ADi7bSPiXr/SOmaolsV/l85j2ciBqiTJkpITE/Dk1nU8vnEVs7t3wKxu7fDX+NHYMfsnrV6azPTYrJ/yLQ4sXYhLe3fi2uEDOLg8VHXcFJ+MnwoAqNqiNQCoqrQqdZ74A7x8/bQG+rY0UsQNMO11IBHlTBxz4qQUCgXkaalwcXWz6XPvnj2JQ2FLUKJ6TdWx6Ht3cHHPDlRr0Qau7u6Qy9MglbrgsZ5XFac3rxU9vmbSWLzfewCkLq7wNFB3AgA2TZ+EdqOMf7jev3DW4Pl/Fs9H5/FTEX3vDlZP1J3mqinjoNjMSklOwv1zp7WOHVsTptPu2YO7ePbgLuJevdI5J0ahUGDDNPFXWEr7Fs83eF5Myeq1AACtBo1Aq0EjAGgvCliiWg0MXbJaK9koUzfY7OcQESkxOXFSW2dOxZOb1zHo9xUmDZJMSUrEosFCZc3hKzZk6rfW6Ht3sGOW8Drm6sF/tM4dCvsDh8L+gEeuXEiKj0eHsRPwT+hcs59xZNUyk9vumjvDaJttP08zeP5h+AUAMJqYWFLi27fYPf8Xk9uLVUAV8/T2DbN7jxLj4wyer9S4mcHzHcdNAqDbCyKRSBBUuRpeP32sVWmXiMgUfK3jpP67dB4piQm4f/6MSe1vnzqO5IQEJCckmPXBqMmUD/Ck+HgAUCUxziDjlGNrM1QqXYybh/bYmIt/b1d94EfeuYVHV8MRF/Ma67//xuxYnt2/p/dcQKkyaNZvsOi5Uau3YeSqzShTp77e6z+d9CO+WLDc7JiIiNhz4uRSkpJMapeWov7t+/bJo6jbsYtZUzdjnz8zOzZnsXTEFzZ9nmbFN2sAACAASURBVL4xOvpE3r6htX/ozyW4dvgAes+cj7WTx2Uplk3TJ+k912HMBHj6iPfKubi6wsXIPx8SqRSuGgXSPH18MxckEeU47DlxQprluk0t6pX0Ll5rf+M006uUKhQKLBs50OT2ZNiF3duyfI8Xj/7D7O7tLRCNrmZ9B1p0JWUXVyGJKVS6rEXuR0TZH5MTJ5SarO4t0Zymq/Qw/AIi79zSOnYqwyBU5ViDxzeu4tl/93HnzAksGzkQzx6ou/kVcjkOLAvFsbV/QqEwXMskO3uv5YeZvja4y2cWjMT6+s0ORc22HVG/c3eL3bP9mPEoVLosmvaxbQ8VETkvvtZxQpo9Jxm9ffUSm0O+ByAsoqaUKvL6Z/HQvjpTabf9PA1fLFgGF1c37J7/C26fOmahqHX1njkfq74ZabX7W0rTvoNwef8es65x9/LCBwO+RHJigpWisrxSNetYZZHF0rXqoXQt80rLE1HOxp4TJ2RoQOXbV9prxsREPYXsxXPRthkTEwCIi3mNuT074e3rl1ZNTAAgX6D6gzB/UHGrPsvc0uYAUKVZSwxfsQGubuZP1x4RthEVGzfTGczqyNqPyb4LEhKRc2HPiRM6t2Oz3nMSqXpKZ/T9u1j3/deZWrsm1kg5eUtwcXXDoN/DoFAo4OXnh/m9O1vtWYXLlsPTWzeMN9TQeshXJrUrVrU6Iq4KReEqNHwfjXv0VZ0z1MtlbYXKlEP0vTsmt89MEkZEZA1MTpxMakoKzmaosCp7+Rx++YXBi5ofwKsnCOuhaM7UMZWp1UOzyjdffps8x1wuZnxQR1wNxyfjp0L24rnO+BRDU22V8gUWs8pqxe1Hj8ejK5cyVXiNiMie+FrHyYiNHdGcTmxOEbOcRCKx7P/qpWrVVW17+fqhZPVaogNnTVk7pvPEaajf+TMMCg2zZIjw9s+Nyu+3QLEq71n0vkRE1sbkxMmkJCfqHHNJH09hzto0WfXVqi0YumS1Re9ZuEx50eOaA3tNNeLPjRi1Wj1lV/byBVoOGq7abzHgS4PXS6C/gu7QP/5CB43xGbXbf2LwXp2+/d5gjQ/fvPnRsGtP+Oa1TC9ScJfPMHDBcri6uUHq4oJPJ0/HiLANGL1mO0b8Kb6uUfXW7USPExHZA5MTJ5MYp1tu/O2rFwDUr3Gsrdv3M+Dq7p6lJe7FdBo/ReeYcrE5pXLBjdH1O+PVZ909vVT1NQBA9uIZgipVVe2XqlHbYN0NBbTHilRv/REAoO3wsfD2zw2pRnExY4vxlapRB8OWrUW3qTNRpFxFo7FnVYNPe+rUKHH38obUxQXunl6i19Rq29HqcRERmYrJiZMRGz+yN3SuaNJiLYGVqqi2lR/wgxetFG1bo017BJQqg0JlymkdL1mjtk5bL5HehbxFAgEAvWbMQ92OXdB68AgEVa6m1WbU6q0YEbYBrm7uAPTPzNEcRyKRSvHZj7/iw2Hqkvwurq7oPH4qfPLlR5dJP2pd26L/UIxZtxMVlWvNaKwlExMdKfq8jAIrVEbknZsG23w6+Sc0+qwvRq3eKnreElVWBy5cobUvdXFF7kKFs3xfIiJL4YBYJ3PrxBGdY7IXz7FwgOWKZhnSf94fWvs9f5oDQJiV4u2fW6cXofnng4UZKwoFZn/WQXW8zZfivTyDF63EyY2rVQsL5ilcFAAQULI0AkqWVrUrF9wYd04dg0+evHBxdYOLqxu+XLoGEdevoFhV8TEWWis4SwCp1AWVmjRHrtx5sX/pArT8YjiKV6uOwb+HiV6vubid5rbCzLVylPIWDdIpdlasSjUUq1JNzxXAoIUrsG7KN3j+332dc5WbfoCAUmWMPtcvfwGt/YqNm5oWMBGRjTA5cTJZLX3um78A3r58kenr8xQqInpcIpFgcOifWPn1CJ2ZJxKJRKunAYDeV0I+efKiQoP3VclJifdqiLZrNWgEChQrgZLVa6mOuXl6orTGQFVN7l5ecPNQT6mWSl1U28WrVccX85eKXmeK3IXF/5uI0ZyZ8/nsULOe0/zzwXDz9ETvGfNw+9RxnVWZ2wwdZdb9lJr15dIERORY+FqHLEbq4oJSNeuo9jO+uuk+9WcAQGszPkQ1kwhNHt7eqP9JN6M9Be5e3gCAwIpV4O7ljUbd+yC4y2fIlTuPyTEY4+7pbXLbxj36ARDqoZjjq7+2okYb9Vo6mmNeMqN8gyaqbeV/IyIiR2HWv3ChoaGoVq0a/Pz84Ofnh+DgYOzZY7is9+bNm1GpUiV4eHigUqVK2LpV/F062Ya+GTFKBYqXVG3nyp3H6KyWjOp16qba/u/Sea1zRStUwtj1u1Cl6QdG7qIejCrJ4odwz5/moE7HLqqEqF6nrmjwac8s3VOX6YXWSteqi8GLVqLt8LFmPUGnQJpUuydKc/CvKZQDfAHtV1RERI7ArH/5AwMDMWPGDJw/fx7nz59H8+bN0bFjR1y/fl20/alTp9CtWzf07t0bly9fRu/evdG1a1ecOXPGIsGT+T4YOAy+GcYcaFL+Ri+RSjF40UpUb9XWrPt7eGf9t/C8RYMgkUot8ht93iJF0aRHP4vPLNKkkJu3KKJPnrwmJV2G6pNo9igNXbIaX/1lXtJvbjJDRGRLZv0L1b699hLt06dPR2hoKE6fPo3KlSvrtJ87dy5atmyJ8ePHAwDGjx+PI0eOYO7cuVi7dq1Oe7KcIuUq6swMadrnC3j5+OqszdNx3CTcOPYvarX9GIXKlIVv/gIIqlgly79Rtxw43HgjET558qLvLwvh7i0+7dVRlKheCw/DL6hn8FhYiwFDsSXke9T9+FOdc0XKC1OSvXz9rJp4ERHZQ6Z/fUpLS8PGjRsRHx+P4OBg0TanTp3C6NHaszJat26NuXPnGrx3UlISkjSqnspkssyGmWN99sMvAIBZ3dTFtQqWKAUA8Mzlg/iY16rjZerU1yqzXlHPeIgi5SuZ9OzRa7YjJjoySyvcai4K6Kg6jJ2A+JgY5A4oZJX75y0SiC9+E6/46+Xji2HL1sHV3T1T9zZUZI6IyN7MfqF/9epV+Pj4wMPDA0OGDMHWrVtRqZL4h1Z0dDQCAgK0jgUEBCA62vCiciEhIfD391d9BQU5/geVrRSrWl21XalJc9E2+sZUuKTXAfnk2ymqY71CDCeKgFB4LLBiFZMXwpO6uGQpMXEWbu4eVktMTOHp45Pp5MSvYIDxRkREdmJ2z0n58uURHh6ON2/eYPPmzejbty+OHDmiN0HJ+GpAoVAYfV0wfvx4jBmjLo4lk8mYoCgp1OMblEXHMgru8pno8cLphdD8ChTE6LXb9c6Eyahi42ZWe3VB9uHt54++vy6Em4envUMhItJhdnLi7u6OMmWE6Zu1a9fGuXPnMG/ePCxevFinbaFChXR6SZ4/f67Tm5KRh4cHPDRqUpCaXGPwZYOuPXHl4F7V/ofDxxpc5VdzEKapiQllX/mDits7BCIiUVkesq9QKLTGh2gKDg7G/v37tcad7Nu3Dw0aNMjqY3Okd7JYPLlxDQDQduT/kCt3HvSbHYqXEY9Qrn5DTgklIqJswazkZMKECfjwww8RFBSEt2/fYt26dTh8+DD27hV+e+/Tpw+KFi2KkJAQAMBXX32FJk2aYObMmejYsSO2b9+OAwcO4Pjx45b/SXKA0IHqsSSu6aXY8xUNMji+Q2zWDhERkSMza0Dss2fP0Lt3b5QvXx4tWrTAmTNnsHfvXrRs2RIAEBERgaioKFX7Bg0aYN26dVixYgWqVauGsLAwrF+/HvXq1bPsT5ETSU3rJWk1ZCQ8vHOhYddeVg6IiIjIMiQKhcL08pZ2IpPJ4O/vj9jYWPj5+dk7HLtQKBSY3V1dZ+bjryejdC3Tkjy5PI1jTIiIyOYy+/nNtXWchDwtVWs/LSXF5GuZmBARkTNhcuIk/gu/qLWfomcQMhERkbNjcuIkIm/f0No/tWmNnSIhIiKyLiYnTkJZel6pVK26doqEiIjIupicOAk3T+1KnsFdetgpEiIiIuticuIklCsJ++YrgAHzl8LLx9fOEREREVkHkxMncePoIQCA1NXFrovNERERWRuTEydx79wpAEDsM8MrOhMRETk7JidOQKGx2B8REVF2x+TECTy5ec3eIRAREdkMkxMnYE41WCIiImfH5MQJOPziR0RERBbE5MQJaI45+eyHX+0YCRERkfUxOXECifFxAICAUmVQpFwFO0dDRERkXUxOnMCeBbMAAM8e3LNzJERERNbH5MTBvZPF2jsEIiIim2Jy4uAOLgu1dwhEREQ2xeTEwT1/eN/eIRAREdkUkxMH9yY6yt4hEBER2RSTEyIiInIoTE6IiIjIoTA5cSJ5ChexdwhERERWx+TEiXw08mt7h0BERGR1TE6ciJefn71DICIisjomJw7Oy1edkHj5+dsxEiIiIttgcuLgPNOTkzJ1guHm7mHnaIiIiKyPyYmDi495BQCo3f4TO0dCRERkG0xOHFhKYiKSExIAABKJnYMhIiKyESYnDuzqoX2q7bSUFDtGQkREZDtMThxYgsaKxPI0uR0jISIish0mJw7M08dXY9vHjpEQERHZDpMTB+bi5q7aLlCipB0jISIish0mJw7sZcR/qm2p1MWOkRAREdkOkxMH5uGdCwDgqtGDQkRElN0xOXFg72QyAEDpOvXtHAkREZHtMDlxYNfSpxLfPnnUzpEQERHZDpMTIiIicihMTpxAo+597B0CERGRzTA5cWA+efMBAIqUq2DnSIiIiGyHyYkDi3stLPrn7Z/HzpEQERHZDpMTB5WqsZaOQp5mx0iIiIhsi8mJg0pNSlJtK+wYBxERka0xOXFUEvVmnsJF7RcHERGRjZmVnISEhKBOnTrw9fVFwYIF8fHHH+P27dsGrwkLC4NEItH5SkxMzFLg2Z3sxXPVtosLS9cTEVHOYVZycuTIEQwbNgynT5/G/v37kZqailatWiE+Pt7gdX5+foiKitL68vT0zFLg2d2qb0aqtiVSdnAREVHO4WpO471792rtr1ixAgULFsSFCxfQpEkTvddJJBIUKlQocxHmQGmpKcYbERERZVNZ+pU8NjYWAJA3b16D7eLi4lC8eHEEBgaiXbt2uHTpksH2SUlJkMlkWl85SdQdw6/KiIiIsrNMJycKhQJjxoxBo0aNUKVKFb3tKlSogLCwMOzYsQNr166Fp6cnGjZsiLt37+q9JiQkBP7+/qqvoKCgzIbplFw9POwdAhERkd1IFApFpmaqDhs2DLt378bx48cRGBho8nVyuRw1a9ZEkyZNMH/+fNE2SUlJSNKYSiuTyRAUFITY2Fj4+fllJlynsnPuTNw5dQwAULJGbXzy7RT7BkRERJQJMpkM/v7+Zn9+mzXmRGnEiBHYsWMHjh49alZiAgBSqRR16tQx2HPi4eEBjxzce6BMTAAwMSEiohzHrNc6CoUCw4cPx5YtW/Dvv/+iZMmSZj9QoVAgPDwchQsXNvtaIiIiyv7M6jkZNmwY1qxZg+3bt8PX1xfR0dEAAH9/f3h5eQEA+vTpg6JFiyIkJAQAMHXqVNSvXx9ly5aFTCbD/PnzER4ejoULF1r4R8keXkc+tXcIREREdmVWchIaGgoAaNq0qdbxFStWoF+/fgCAiIgISDXqcrx58waDBg1CdHQ0/P39UaNGDRw9ehR169bNWuTZkDwtDeu+/9reYRAREdlVpgfE2lJmB9Q4m4t/b8ehP5doHRu7fpedoiEiIsqazH5+s/SoA8mYmBQoUcpOkRAREdkPkxMH9tGIcfYOgYiIyOaYnDioTydPR77AYvYOg4iIyOaYnDioYlXes3cIREREdsHkhIiIiBxKpirEkmXFvX6Fszs22TsMIiIih8DkxAEsHtrX3iEQERE5DL7WISIiIofC5MQBFa9Ww94hEBER2Q2TEwfUtM8X9g6BiIjIbpicOCAv3+xbop+IiMgYJicOSOrKccpERJRzMTkhIiIih8LkxM7unT+jc8zDy9sOkRARETkGJid2tv2XH3SOSV1c7BAJERGRY2ByQkRERA6FyQkRERE5FCYndpSanGzvEIiIiBwOkxM7evv6pc6xdqO+sUMkREQ5mFwOTPEXvtJS7B0NgcmJXb149J9qu0qzVhi9ZjvKBze2Y0RERDnQ32PV24dD7BcHqTA5saOds9V/CVoOHMZZOkREtnZzJ3B+uXr/2CxgemEgOd5+MRGTE0fBxISIyMZinwLre+keT3kHXAgz7R7yNODZdSBe9zU9ZR7rpBMRUc5zYApwfI7+8/9MAMq2BvKX0d8mJRGYHqDenxJrsfByOvac2JGrm7u9QyAiynl2jDCcmCgtqAWcDtV//tFx7f2Xd7MWF6kwObGj8g2aAAAKliht50iIiHKQiyt1j334i3jbvd8CSW8BhUL7eKIM+Kuz9rGVHS0THzE5saeYqKcAgEpNmtk5EiKiHK7eIKDdXMCnkO65kEBgam4gLVXYP/U7MCNIt53sqW4S4wyubgJO/ga8uG3vSFSYnNiRQi4HAKSmcF49EZFJ3r0Gzq8AEt5Y5n4VOwAf/ixs1/4cGGfgA3pJU2EA7D/j9bc5PkddM+XVfSAlwTJxWtPFP4F9k4Doq/aORIXJiR25eXoCALx8/ewcCRGRk1j1MbBrFDCzeOau3/+devuj2UC3VUC9wdptvn8DdF2le230VWBNV93jvkXU2wenqrd/qwlMLyQkNICQqCi3HcXtvcB/R4VtD1/7xqKByYkdKXtO3L287BwJEZGTiLqs3n4TYf71N7artyt3Em8jkQCVOghJSkb3Duge+/Kk4WdOywtsHSokKtPymh6rpaVmWDLl/HJgbTf1vnsu28ZjAJMTO0pJTgIASKX8YyCibCo1SZ1QpCYDJ+YD55YBSXHCwNS456bfa11P7f25Vc2PJ+ahetvbSKIgkRg+P/QkMPkV4JXH+HMvr1Fvm1Ii/+FxIOaR8Xam2vMN8GMBYOdXwn7yO2DXaO02AZUt97wsYp0TOwjf9zdOb16L+DcxAAAJkxMiyq5+LCh+/OgvwNsoYXvyK8DFhI+jW7t0jyW8AbxyZz4+Yz4OFQbAPsswHqPnZu0P8x4bgTWfmnbPC2FA3YHi565sALZonBt9HfAPFLbvHQRkkUDN3iaHr3JmkfrZ7ebq9gBV/dS0JMtG+KloBweX/a5KTABAImV1WCLKhuJf6T+nTEwA4Id8xu/19KL48S16PuQzenJBGKSq1GySaddV7wEMPQ4UytBLU/YD7f1yrUy7HwD8PU7oUUp/ta+SKNP9eSJOC9+T4oC/PgF2DBd+jjv7hAG3plj1ifb+1NzAhgwJTvv5psdvA0xObOz2qeM6xyTGug6JiJzRL6VMbxu+xvD5JXpKLtzdp54do5xtcnsPsKix+sM7JRFY2lz7uuINTI8NAIYcB8beAXpuAibpeRXVWmPRwCmxwHev9d/vx4LANI2eipRE8enJmwcI6/w8OKR9fM2nwoBbU9w/aLyNm2ONfWRyYmO75s7QORb7PNoOkRARWZEp4yo0bRtqetsJUeLHFzUSCqOt7Q5EXxE+vE8t1C4xr4ovWfeYMb4BQNmWgKuH+PngL4FvHgLfpfeMS12AltMM33OKv9A7Ihaj0sMT4msAAcLYEblctxdGaWF9w88HgPFPjI+vsTGOObEh2UvxbJtjTogo29GcsptVqUna++7e+ttmHEvxzwTdNh7+5vecmCrjuI36w4B8ZYGgesDtv4XXMhktb234nobGsvxUWL3deytQWqOH6OEJ4MVN4zE70BRiJX4q2tDOObq9JgAggWNlrEREWZL8Djj9u2ltOy5UbysrsGZ0bql6+8v0MRhNvs5cbADw9X39vR+W5uIKVGgL5MoHxJsxM6mvyOBfY1Z1EqZXK19zhbXVPj/sHJCrgPaxTn+Y/xwbYHJiQ8//eyB63N3bwG8BRETOJOaR9m/zADD0lHq7wUjtc4F11dt7vwGWtRYKgynLwMvl2r0ffkWF781EekSMafsrMO4u4OJm/rWWUHuAdixiSjUDRlwESjYWP/9djPZ9MtI3vfq7GKBAOWCMRk/K2DvAe93E29sZX+vYyIXd2yDX81tB2XpW6l4kIrK1zRk+OOsOAgIqAaOuAi7ugG8h4KTGzBDlNFlA3UOy1sAHpmd6RW2JBBgZDlzfInxYm1IxVt/0XVvxyi0MlAWEMvx/j9Nt0yYEyJe+GOywc8DCOtrnpVKg3WwgfzkhmTNFUD3hOkBIzJQxODD2nNjI4ZVLRY+3GjISbu426l4kIrK0FW2FVwgrPwZmlgQKV9c+3yi90FfuYkJikpGh8SMZZXwFkbck0His8KHfR6Pya+dlwgewZll5R+OppzaLv8aMnQLl1INrASCvxgr29YcIP6Pyy5AB+zIfp52w58TOytQJtncIRORMUhKBK+uAoPpAwQribVKTgDmVAYmLUMTLlAJnmfXohPBdOdX13BL1uQrtAD8LJgiGXkGUaiqUm9ecdTLmhjBr5+LKrI1RsQapVFjb5/6/QlG0l3cAN0/Aw0e33aTnQp2XoLri9wKAAhXFB79WNbEwnINhz4mdefk43ihpInJgB6YIJch/rwckvQVu7hQSFk0/FgTiXwBx0cD+yZZ9flKcejzIwxP62/kXA7qvFj83IRLovhaYmF5GoZiFXm1nnA4rkQCF3wM+miVMA3Y0dQYI/418CgAlGgJFa4m3c/UAigcLU5P1GXpCWGG56qfA6BvCsaD6QGfxXntHx54TO2o5aIS9QyAiZ5IQA5wJVe+HaIzX0BzLoOn078I4Bk0v7gi/oWv2asgiATdv/aXg370Wnv9bTSGZ8A0Arm81EKuBAmTuuYQZLEq9txqu80HGSV2EFZaVnGBciSFMTuygaIXK6PrdT5C6sGw9ERnx7jWwfZhQI8OQJc2Bgf8CV9YbbvfqvnqQ5ZATgG9hIDURmFNJODbuLuCTYT2cy+uBrYPU+xFGVuEFgNr9jbdRcvMEPvwZ2PO1MF7k0l9Cj0JaMjCzhNAmqJ7p9yOnx9c6dlCkfEUmJkTZ0dklwuDQa5t1z/39P2CNxpiJUwuBW38Lr0m2DgFe3gVePxDWV9H0c0njiQkAPE1fO2bvt7rnNKu1Pj6r3l7UUCgxf2ev+tivZYUE5sAU4btCoZ2Y6NNri/Z+7mLGr9FUb7Dw237VLkCfbULvimZBs64rzbsfOTWzek5CQkKwZcsW3Lp1C15eXmjQoAFmzpyJ8uXLG7xu8+bNmDx5Mu7fv4/SpUtj+vTp6NSpU5YCd2YPLpxFkx797B0GEVlKajJwc4d6auim/kC+MsIKsE2+Fha5O5s+0+TUQiEJUU6bdXEXeggur1XfTzmwUzm2I6v2TVa/DqrYQff87jHa+8o1W47PMf0ZpZqqt6t1N6/nxJDv3wjJlau7Ze5HTsGs5OTIkSMYNmwY6tSpg9TUVEycOBGtWrXCjRs3kCtXLtFrTp06hW7duuGHH35Ap06dsHXrVnTt2hXHjx9HvXo5o5suMT5Oa79iYz0LWBGRczocAhyfrX1scRPh+/nl2sczllMXW+Nlam5gzC3gspHF8AAAEgBGkhjNcSo3d5hwTzMNOiKMeVAOcLXkInISCROTHEiiUGQ+NX/x4gUKFiyII0eOoEmTJqJtunXrBplMhj179qiOtWnTBnny5MHatWtFr8lIJpPB398fsbGx8PPzy2y4dvPfpfPYMmOKar/zhGko8Z6Jq0kSkeOblh+Qm7nQnbnqDxPGk7x7CXz6J7Cxr3D8u9fAtLy67fvsAFaK9JJkVZEawOd7hAQk7rlQDt3BFo0jx5HZz+8sDYiNjRVGA+fNK/IXI92pU6cwevRorWOtW7fG3Llz9V6TlJSEpCT1Qk8ymUxvW2fgkUt73nreooF6WhKR00lLtX5iAgBtfhK+lCrrmY3h5g30/wcoXM3w/cq0FKqzXlhhXhw9N6t7RjIOnCWykEwPiFUoFBgzZgwaNWqEKlWq6G0XHR2NgADtKWIBAQGIjo7We01ISAj8/f1VX0FBQXrbOoNja8K09n3z5rdPIERkefEvtPe99P+ylmma1U/FDL8grEY74iIwMUqdmAQaKNr12TqgUkf1ftPxQEnxHnAAQLFgoOkEYQE7IivLdHIyfPhwXLlyxaRXM5IMXX4KhULnmKbx48cjNjZW9fX48ePMhukQnty8ptouVbMOJFJOkiLKNpI0enar91KXa89oUoYVaf8nvhAoxj/RPVaqqeEY8pcRaoXkK619PGMs1bqnP+OpUDW2dDNg2FlhrEjTb4FuGYqmKUunB9UH+u8Fmpq4lgtRFmXqtc6IESOwY8cOHD16FIGBhl9RFCpUSKeX5Pnz5zq9KZo8PDzg4ZE91pvJOKTn46+/s1MkRGQV9w+pt9v+DLh6qquyFqoG5MovrIni6gFMfiVM2y3eAPDOqy6Udf9fYbn7cXcBD19hhspUPcXQzOGj8e+siwfwyWLhS1MBjdmWnn5AzT5CufdvI4TS6U5ezIuck1nJiUKhwIgRI7B161YcPnwYJUuWNHpNcHAw9u/frzXuZN++fWjQIGesxBtx7bLWvqEeIyJyQporw7qnz1qcECWUlS/zgfZrEBdXoGI73XuUbq6dBGj+O9E6RLe9qYrUEL77BwFfXTbcVqnDb8IXkR2ZlZwMGzYMa9aswfbt2+Hr66vqEfH394eXlzBAqk+fPihatChCQoS/UF999RWaNGmCmTNnomPHjti+fTsOHDiA48ePW/hHcSxpqal4+fgRNv04yd6hEBEg1CI5MkOY9ZKVcRPXtwll3NNShMXaxLh7G16kzhRdlgMv7wF1TSiApo9UKvTCKOSG12UhcjBmJSehocJc+aZNm2odX7FiBfr16wcAiIiIgFRjTEWDBg2wbt06TJo0CZMnT0bp0qWxfv36bF/jZM+CWbh96pi9wyDKOZLeCqvxytPEF3n7sYDw/dgs4fvo68JsI3YTKwAAIABJREFUFXNc+BPYOVL8XFB98+5lTJXOlrmPRCKsTkzkRMx+rWPM4cOHdY516dIFXbp0MedRTo+JCZENPb0grC2jqflkoGAlYYG5rUN1r5lTWf0qJeI0sLy1sHptnS902yoUQGIscGmV7jmlxmMzHz8RaeHCfzbk4uZm7xCIsqeMiQkA/PuD8L1AReDFTfHrpvgLC81tHiDs7x4LnA4FvjwjjA8BhFLz82sYj6FsS/PjJiJRnNNqBWmpqaLH244YZ+NIiEhvYqKkTEyUXt0DfsgHyKKEfVMSkz47WCWVyILYc2IV4q+/ytVraOM4iLIxhQK4dxB4fd/0az5ZApyYBzy7Zrzt7ApALgMVUHttBvyLAQXKmf58IjIJkxMick7/HQFWZxg0OuoaMK+aMDtFTNVPgXevgL3fmvaM+Ofix5tNEqYJE5FV8LUOETmn9X10j3nnA0aGi7evN1R49VJ/qJCkZNRwFJC7mGnPfv9/psdJRGZjzwkROackkcql7t6Am0aC0XEh4FdU6C3RnJrbeSlwdaN6v9yHQMupwtfNncD6Xrr3HnJcmNVThCuKE1kbkxMrSE3WXaHU09f0paKJSERCDHB+OVDpY8CviO75ISeE7xKJsCrvy7tADZEkQ8m3MPA2fdBrp0Xq4xXbC1OMp/hrtw+oAhSqmrWfgYhMwuTECjQX+lPq8/N8O0RClE0kxQEzSwjbB6fpns+4/kux+sKXIUNPAj+XBCp3ArxE1rEZehIITV9mY9JzzsYhsiEmJ1aguX5OkfKV0H3qTK6pQ6SkUBj/oJfLgdREoeT6b7WB2Aj9bav3zFwcmgvviQmoDAw+CkjdhEX7iMhmmJxYg8a/ux//bxITEyJASEqUK+0OPSl8+ItJSwF+yG/6fdvNyXps+hR+z3r3JiK9OFvHCqQu6pzPlJL/RNmePE27imtoAyDyknjbk2a8Av0uhr0aRNkQkxMrUKSlqba9/fwNtCTKIWIeApEXtY/90VS8baSeqcCA9kyZya+EVXeJKNvhax0rUOipEEuUY93dL35cOSOm2URhteDURPF2PTYA5VoL21FXgDzF1WvfEFG2w7/dViBPFXpOCpVhWWvKwe4fAlZ9DNQdDJxdbLjtoenixzsvA6pmWNG8cDXLxEdEDot9olYgTxMW/nNx5SrElEPJIoXEBNBNTEwtYlarn25iQkQ5ApMTK1CuSuziyo4pyqFmVxQ/3ngcMOiQafdoM8Ny8RCRU2FyYgXy9AGxTE4ox0lNAm7v0X++yTjh+wdTDN+n9zbAzctSURGRk+GnpxU8uHQeABATHWnnSIhs7MeC+s+Nvq5OOBp8BRyYImz33ASUbQkkxwtr4Ji6+B4RZVtMTqzgzqljAIA30VF2joTIhuRy8eOTnuvWIpFKga//A55dA0o0Fo655xK+iCjH42sdIsq8Pd8K04Ff3Qce/Kt7vu4g/UXSvPMCJZtwzRoi0sGeEyIyX/RVYOXHwLuXwv5vGWbgTIhMX5PG3faxEZHTY3JiRb75Ctg7BCLrWNTI8Hm+niGiLOBrHSsoWaM2AKB2u4/tHAmRBSXKgFWfAJf+Mtzu41DbxENE2RZ7TqxA6uICAHD14IJklA1EX9XuKbl/0HD7ql2tGw8RZXvsObECZZ0TqdTFzpEQZZFcbvgVTusQ7f2K7bnmDRFlGf8VsYKou7cBAFIWYSNnt7Cu4fPBXwKPTwM3tgODj3HdGyKyCH56Wtib6Cgkxr0FACj01X0gchav7uo/1yp9sb6uK20TCxHlGHytY2HLvhqo2vbNl9+OkRBlQcIbYEEd7WN5SwH1vxS26w4CGgy3fVxElCOw58SKXjx6iGJV3rN3GESmO78c8PADNg/QPt7/H6BYfWG7TYjudUREFsTkxIokUnZMkRO5tRvYNVr8nDIxISKyAX56WpFySjGRw7u2GVjXQ/xci+9tGwsR5XjsObGQF4/+w92zJ7WOefr42CkaIjOs+Ah4dFz83MBDQNGa4ueIiKyEyYmF7J7/C149idA6VrqWkWmYRNZydROQkgDU7G24XUqC/sQEYGJCRHbB5MRCMiYmAODi5maHSMhiYh4B89LrdkyJtW8s5ghrBzw8JmwXrQkEVNbf9qWBqcLfPLJsXEREJuKYEyvK9hVi51YFpvgDT87bOxLrWN5avR3/SrxNSgKQlmKbeEyhmZgAQGgDQKHQ337nV9r7XZYD458A3zwEvHJbJUQiImOYnFiAslx9tqdQAHf2AXEvgKcXgDfpvUVLWwApicIHdXbx+gHwNkq9f3y2bpukOCAkCPipqG4CIE8TErcp/kDyO/OerS+ZMJRk3N4DnA7VTkyUpuYW4hETeVG93W4uUKUz4OELeOUxPV4iIgvja50sOrdzC47+tdzeYdhG+Gpg+zDxc9MDhO/lPwIkEqDTIuFDzhIUCuDEPCCwDlCioWXuacgUf91jpxYApZoKYzmurNM9f+RnoOk36v2NfdXbPxUGJr8ybc2ZRBkwI0jY/i4GUE5Hn10ZkD0Bvn0MePppX/P0ArC2u+H7TssLlGgMlG0F7J8MuLgD/kHabWp/bjw+IiIbYM9JFuWYxGSKv/7ERNPt3cCtXUBIoLCvUABnFguvgI78bPi3f30OTgUOfA+EtTX/WiMi4yIR/jzctMaru4gnJgBw+Cdgy2DgXvqKvbJI7fPza5j2jP2T1dvT8gCv7gOHQoTEBNBNQtJSgSXNde/z4S+6xx4eU98/LRl4fV99rvYA3fZERHbC5CQT7l84i1dPH9s7DNuQy4ENfY23E/PPROGVwp6vhVdAh6YL+0/OC8nOgalCu7QUIOI0cHOn7j1e3QeOz1HvK8d+PDqlHvNyc6f+1xYGnIs+h9abW6P3nt748fSPwsFlrQ1fZMiVdcBfnwC7xwm9GZpidQdMi8o4QPW3msCRGer9Rye0z2d8DgA0HAXUGyT0vJiq7a+mtyWibOX4hrtY/r9jSHibbO9QVJicmOnJrevY9vM0hI0ZarBdj+mzbBSRla3sANzYlrlrTy0QP760hfD9+GwgIQb4Ib8w+HR9L+DKBnW7DX2ED2dNv5QSxlasaKMe87K+l/DaQl9isf87IYlJigOe31SNBen/T39Vk/W31wsJzuPT6utKvg98/8bMHxrAuSXmX6OUMfkQk5okfL+9F1jeSvtc0wlAy/SkTyoFCpu4fAKrGRPlGClJabh34TkS41OQ8DYZl/99jIS3KTiy9ra9Q1PhmBMzhe/dpdpOS9U/S6NwmfK2CMe6Ds/UHWD5yVKg2qfAu9dAkgyYl/7hl68M8Oqe+c+YWUJ7/+lFoFpXYfvGdvFr9n4rfvzxafV4EeXU39QkYbwKAIQUBQAoAMRJJDqX754ViMZSCfzk6a+e+u4Qvg/YDyxrqW5YdzBQsR3wZ3v9P5eYxU2AAQcAV3edU4mpiThw7Ee0AqB7NoPL64CafYC13bSP994GlG6mfazvLvUYFn1GXTP2RCLKJuRyBf746ojouVI1Ctg4Gv0kCkVmBgHYlkwmg7+/P2JjY+Hn52f8Aiua1a2dSe3Grt9lvJEjeX4TWPw+0G834OEDFKggvILRFDwcaD1d91q5XPjN+8BU8Vkt5qo/DGgyDvi5ZNbu02sL4FcU+L2e6tBrqRTvFw80eNnV/yKAfn9rD75Nfgc8OScMynX3Fo4lvQVSk4H4F1rP0DLionbvT5GaUAz8FxLN5Oj5LTTa1RmxLlK4KhS49PCx8Jrl73HqNs0mCq/F9Bl8VH8viUIhDFJOjBWSNZ+C/2fvvONjvv8H/rzsvYkgkRBbEntH7L1X1SytUlqqNUppffErrSpVSlW1du2i9ooZasbeYiSyZM9bn98fJ3e53MglEqL9PB8PD/m8910u93l9XlP1O0uLBqeyRt8LERGRfxfndj3kwt4IvX1jl+vxX3tFCnv/FjUnxUDXT6fkP6ikkJUCGfHw88vCbr+11T9u4Bao0l5/X45JoPUMcK2gigqxc4dvK2jGTHuuiloxhbNLVf+09rAApdy0+Tms6w2eAerLffZ2TC7tYdrcvFFBVnZQMUS7zdoRrAF7d5WmRpoBm4fA/cOaMe6VtKYEWMfDmkDqlq7LqkwbLt7dwU0rK5LdVaG78hyhpeFIlXZkTmlw9YMWk4wLJ8bMNzlr2uSKQjIzEwUTEZH/GHFPUw0KJiWNAhuaT5w4Qbdu3ShbtiwSiYS//jLujxAaGopEItH5d/v27UIfuiQzZuUGqjYJftPH0EWaDqs6wbJmKtPH7T2q/+d55x9JUmeIYcEkN2ZmUO891Q3Z1gU+uw1l68K4K6qbu3UuqfnLaJgaCfWGq8JbP79reN2vEuGrFzDmrHb7aP3+GXLgoaVK7j6Sep+7lpbIQK9g8k5Kqu4Co42kczeGlR0M3gbT46DHz6rXDxD0LgBferiph16KvUTt1DO87+XJAnftnCK7u85W/WBhrRJ6xl9RCRi9DfiydPmX+DeJiIgUK5v/77zBvlE/hRjsexMUWHOSnp5OUFAQw4cPp0+fPibPu3PnjpZKp1SpkmPbKirKVKqMreObNTsZ5Js8T8mGKtDqo+Pcwu3p5AUfHtNcD94Of38K7eeApa2qrduifPaep9HMlK6u3Ve6OkyPhYhT4BdCxubBNJLe0FnCQhA0Gok8THuRyKSERGZ4uLPPwR6Azcm36F8mQO/43GTKM7E2t8ZMopHxb7y4wYC/VeG+3dNvMjpwNE8bDGav9DG7sp7luybAtBu/0q3+ON2OwP6wfaT6MksiwWbSQ7Bz0x0rIiLy1pGRIuXExjsEtfHGy7/wGZqz0mXcOBlJQlQ6rYdWx9zCuB6iOMw5r0qBhZNOnTrRqVOnAm9UunRpXFze7nTYcplhB9gek2ZQIbD2azxNAdCXVMxURhwsumRq3g3gIyPRKH1XwVZNBA3lG0Jj7agoAehc3ouqUhmLzMw5H3uJEac/pd79eriUKQNPdIUTfYJJu1L1+faf7ZgB1gJMf5GgFk5mn52Nu607bXzaGDzq6EOjOR2lei2nBpxi3NFx9K3Sl2mnpqnH7Hqwi10Pdhl+vUZourEp2fJspEopV4ZcwTxPKYRhXqW5ZGMDW0K4NuxaofYQERF58zx/kMz2+dopAR5cjmPkohZY2VggCIK2j1o+xD1JZfM3Gg3J3X9i9I5r0MWX83siGDKnSeEOXsy8kkOsRCJhx44d9OzZ0+CY0NBQWrVqha+vL1lZWdSoUYPp06fTqlUrg3Oys7PJzs5WX6ekpODt7f3GHWL//vE77pw5obfvtTnApkTBkzCIugwXV8PEe2Bpo+mPugIrXqrnxl+FfVPg7j7T1i5bF0Ye1fgoFAEF+cOSKWUkR1/FI+IM1BuOYGVPtiIbGwsbErISuJt4l7i4m0y7rMp7Us2tGrcTCm4eHFN7DP2r9Mf924pa7QF+PlrX/ar0Y3KDydhY2Gi1K5QKaq/VCKLW5tZkK7IxlQ5p6Rx4KQgB1FZaEFC1F+UtHJh763ed8e0qtOOHlipH44yMFwze1Yd72ZpaP8Hlgvm57c8m7y8iIvJ6uXc+hoO/3aBKI098AzyoXF+VUTs7Q8bKz/SUnMhDzeCytBxUzaS9lo4+mu+Y16kpKbEOsV5eXqxYsYJ69eqRnZ3N2rVradOmDaGhobRo0ULvnLlz5/K///2vuI9WYAwJJn2/nPP6DrG2F8TluiEfmQUdv1H9HHNDI5iApqKuMSq2hIeh4OIDHxwuUsHky1NfsuvBLmY1nUWvyr3yHV93rSqqZUPnDSTGXWLsEVVG2n5V+rHl7hYAytprzFOFEUwuDb6EpfnLatFtvlZln32Jk5UTKdIU9fWWu1tws3Gjunt1tRYlYLWuuacgggnArPgEvo17Qay5OV4KhVbF4/0vrnA59rLW+JzrdFk6jbe01FnvZORJVl1fxcKLKqHt8pDLWJiJvu4iIiWFg7+pNLp3z8Vw91wMPjXdsba14NJB05Iz3jgZZbJw8m+h2DUn+ujWrRsSiYRdu/SrvEuq5sRQGPGEDTsxM39NFYj1mWhywkjzM9/0XwPbRoJPI3j0UtD6OqlIBZIcrsReYci+IVptHX07Mj9EO6367ge7Ke9YHlsLW/rt7ldk+y9ru4xGXo3UAg/AgT4HKOuQx/fmxQOIvgbVu5Emz6TJxuJXcV7LdldFzrj6QbNx4KbR4MRlxNF6i/ZTTSnbUhztf5Rh+4ZxKfZS3uV0+CDgA8bXHa/TvvbmWr47/x0LWy7E38WfUnalsLe017OCiIhIUSHNlPPrBN0H27HLWxO24z6XDpiYPRoYOLMRrmWM/83mpzkZNKsxLqXtTN7zVSmxmhN9NG7cmHXr1hnst7a2xtra+jWeKH8Ucv1hrBUC67w+wcQQv7QAU24yNXqo/gmCKrlaqerFIpgIgqAjmADsj9jPlbgrRKdHF/meuQkfGq52Us0x/TTxaqIrmIAqsuhluK+DlQOnBpyi+Z/NX2n/HB+QZ6nPcLNxIyk7iZ47e5Ipf1m1+cNQg3NL2ZXi49ofs+TKEsrYlyE6PZq4zDi9GpselXqw84FuorqV11YysNpAErISyFJkEegRyNPUp3x3/jsAJoROUI/V588iIiLy6ihkSi4dfMw/ux8ZHCPN0pTdcC5lS3Kc8cruG2aeY9RPIUTeScLL35lHV+K4fPgpzqVsadKzEg5umvtmzRblKF/VlQO/XsfL35neE+u9+ot6jbwR4eTy5ct4eZmY86KE8PR6uN52G4cichY1hVQjN3VZuubnsnUhKs8T9vD9mp8lEvDTb1J7VcYfHc/Rp4Yl96IUTDZ33Yy1uTVeDl40XN8QUGlHckfPrOm0hsuxl6nvWd+kNZ2tC+88/G3wt3SuqClOWN5RlezNztKOsHfDWH1zNd0rdc93nVFBoxgVNEpL65MXXydf5jSfo1c4AbS0L+Prjufw48N6xy26tIjP63+e75lEREQKxvbvLxL7WE+qgpfk1nA06+uPb4AHEdfi8anpzu2w51w2YPL55RPd7K4vnqXx8HKcVluzvv5YWpnjX6/kReKYQoGFk7S0NO7f16Qpf/ToEVeuXMHNzQ0fHx+mTp1KZGQka9asAWDRokX4+vpSs2ZNpFIp69atY9u2bWzbtq3oXsVr4OJe/TeBhj36vr5DXNB1lsxBAC7aWBOQnY31iAOQHgvXt0H9EUUXbZMPiVmJRgWTwvBL21+ITI8kLiOOZeHL1O3rO6+nursmtNhQxIqthS1NyzYt0J5Xh15FrpRTd51+4eDMu2eQKWWEbNL49+zovgN/V3+Da5qbmTOi1giD/fo42u8owZt0c+ZMqDdBvdbvHX5n+IHhAFR0rsjD5Ic643+89KPBPf648YconIiIFAP6BJMR85uzapJuHqWEqHRqt/WhtqfKKb9pb38q1/fUiropKJZWb7dGtMDCyYULF7QibT777DMAhg0bxh9//MHz58958kQj8UmlUiZOnEhkZCS2trbUrFmTPXv20LlzZ521SzIR4RpNhE+tIASlki7jJ2Pv4mpkVhGTqzrtuSGb+ODUJL6wr86gyDuMsk4nzFaVO+SahRU4l4dm2n4HWfIswuPCqetZF0szS53lw6LCsDSzpH4Z07QMuREEgRabdLUxB/scpP02ExK4vWRvr73EZsby3v73qO9Zn6blNILFezXfIzo9moSsBAJLmeDsW0gkEonGafYlMxrPYN2tdazuuBpHK5Wwt7PHTsLjwunh30NLW1NUuNi4cG3YNS2TjrW5NUNrDFVf1y9TX0sw02f+yY/xR8fzab1PcbJywt3W/dUOLSIiwj+7tR8SrO0tGP5dc8zM9JvRqzUpo9NWyseRMctaIZFIeHgljvsXY7l3Xn9YcF6qNtJd721DrK1jIrmdYSds3InZ67bTC4K61k2imRktctWHKe9QnmdpmgRfhrQIOTeu6m7V2dxNVf1378O9TDmpm27/+5DvmXh8Im192rKw1cJ8j/fX/b+YcXqG+npZ22U0L6fy3ZAr5ZhLzPm/c/+nqv4LNC/XnFORp3CycmJFuxXU9KiJTClTC01KQVksN/yCEJcRx8KLC/ko6CO8nfIpnleM3Em4w7ij4wguH0yfyn20NEZ5SchK0NLo5KVZ2WYsa7uM9w++z/lo3aeyvb320nlHZwI9AlnfZX2RnF9EG7lUwd3zMdw8FUXMoxRqtShHyMC3u1CooBQQwODN979CenI25/9+xI2TUeq2kIFVqdWinPp6+cehKORKrXk5QogxZFIFK8bpmnRK+Tgiy1aQFJOhbsvJkVISKOz9WxROTCS3cPJGivpd2wrb3kcG1M2TjyMvgR6B6rwXIZtCUAgKLgy+QP11Go1IYKlAevv3ZmbYzHy3vjr0qtE/nFRpKk03aptO9AlIgiDQZksb4jLjODvwrBgpUkxkK7KxkFigRKnjt3K472E87T3Z83APX5w0UN05DwtbLqRtBQM1l0QKhEKhZPnYUJ32knQzKSiCUmDrtxeQSZX0nVIPaaYcB1eb/Cf+Czn8x03unNX2q8sreKQmZLFm2hkAHN1sGPqN6WbnI6tvcjssmmZ9/bGwMtcSekI33OHGiUi6fhxEhVolRwMqCifFzBsXTr6vAmkxfOvmwjrn1/sebOyykVoetdTXByMOsj9iP9+HfM+KqytYemWp0fEib44b8TcYsGeA+jpHaJQr5dRZm09NpTxcHHwRK3OrIj3fv52sNBnXjj/DvawDfrU9+PmjY0bHd/0kiAo1S86NxRRy32xz6DYuCJ8ab9frKAryhvH2mVyPMhV1nexfRKaRFJtBpTqli3T/gmaTfR28VaHEIgVEECAtRieD6dAaQ1lzc02xb//unnfVN7UUaQqfH1c5UAat0V8Jt4Z7jWI/k4hp1PSoybmB5/jh4g8MqzlM3W5hZqHjz5If9dapQhHz06SJaFg1+RSC0vTnv79/0kQFWlqb02lUAHKZAr+gkluLLHS9bjLE3YvDGfVTCBaWxWv+zkyTcuv0c+q080Hyhk1KSj2/Z1cv/dph93IOuJdzKPIz/Jv+Lt+sUf8t4U0qlwRB4JttPXUEE4CJ9ScyuPpgHK0cuTr0Kr+0+6XQ+yxpvYRt3bdx8p2TVHGtotO/LHwZgiDQbGMzg2vs672Pa8OuvXFfERFt7CztmN54Ot6Oun4zM5vMBKCMvekOdIFrAhl9eHRRHe+tIStdRuj626S8UOWiSIrJ4P7FWKPfDwURTPIiy1awa/EV9i67xh9TTr3R7yFjPLmRoLd923cXiXmUovemDSrfm1d9TasmniJsxwN2Lb5i8pzIO4lcPfaUK4ef8Ounx7l99vkrnQEg6l4iy8boasWsbcXn/8IimnVM4NyOzZz6U6WhqNupO63e+/C17W3oyXZqw6kMrK5bWTg+M55Wmw3XLQKVT8rV+Kvq6z+7/klN95om762P71p8Rye/gheEFCkZCIJAr529eJD8gO6VunMp5hJ1Pesyo/EMGqxvoHfOf6ng4Pk9j7SSaQ39pqmWKWP00paYm6uEckEpcONUFMc33DG4nm+AOxHXXhjsN8a7XzfCzcAT+evmzrloDv9+0+gY7xpudB+nXRT15Ka7XD2mcuI3VudFmiVn/ddnyUiWUq9TBRp1r6jWDjy4HMv+X66rx+pzKn0UHsfeZarPaYsBVYiPTONmLmfVHHp+VodyVXQjL2+cjCR0ver3OHh2Ew6uvE7s41TKVnah1+cqfy5ZtoIV47UdVQNbl6dGs7LFoh152xB9ToqR3P4m1ZqF0GXcpNeyb7YiW8uJNTfejt7s7b3X4NwzkWcYdXgUAA3KNFBHZuSOoll5bSXuNu4G694Y2z83e3rtwcfJuJOuyNuLIAgErtEN3f42+FsCPALeaCTT6+Dx9Rf8vUR/EsYcSvs60ePT2vz6qf76W7mp3c6Hxj0rkhCZjoe3AwqZEgsrc72+G4YoiMlEIVPyz98PqdmiHE7utibNyY9bZ54jlyo48edddVuvz+tQtrKr3vTpuQUQQRC0fG8MpVPPK3zkXSvvPkGtvWnev7JWmylF8HLo+nEQfy8J14quyW9+v6n1USoEtn2nXVX4dRbWK+mIPievCd8gw1k7i5LdD3Yz7dQ0g/2tvY1/+JuWa8qhvodwsnLCztKOI0+OEJkaqRZMQFWDxRjW5taUcyhHZFqkTp+bjRt9KvdhXN1x+bwSkbcdiUSi1z8lJwR9WI1hTGww8U0crdhJisnIVzABiI1IMUkwAWjSsyJm5maU8lHly7F4mSzL0c2GYXObYmFlTkaylI2zzhk5VyYe5U17Kl/+SSgAlw48YfRPLTG3fDWzq6Ebtpe/KtWBraMlmakyrb74Z6l4lFe9Xlm2Qqtv/VdnadrbnzPb79Pj09qUr+ZGUkyGXsEEQCFXkhyrm+Y9/OhTLeFEqVDqjDFGzu/5+IY7WlEwxtgy94LWddXGZWj7nuhzVxSIzgEFpEaL4peIk7OTdQSTpdGxDKisKYxnys2gjH0Z7CxVTyRtfNowtObQfGbosq/3Pr3tv7T7RRRM/mNcHXqVjV02Ymuh/fS9+ubqN3Si4mf912dfeY2uHwcxeklLWg2pxgc/BGNmbvhr18HVBht7S9zK2jN2eWv1v7zIZaobvCAIZKXLdPoNsfyT0ALftHOTmpClt710BUe1SeWd6Q11+jfNOU96cjYyqUKvEHdmuyrr+M5FV4h7kmr0fV/+cSi7fzLsY5KRIiU5LpPUhIJVC8/N0tFHibqXVOB5bYYZzkEkUjBE4aSAvA5v6JabW2pdz42Np0VmFl82/Yprw669Vlu/RCLh5Dsn6VKxC7930KTPr+b23yrfLaL6LNTyqMW3wd/q9EkVUgCi0qKYeHwiZ6LOsOPeDhRKhc7Yt4W8N/Gun2hHp3lXd6XLGOOZikcuakGFWu6YW5hRo1lZrO10MzObwtjlrRmzTONLtu3bi2rzyG+fn2RCgcmzAAAgAElEQVTp6KNEP0wmPUlzQ455lKJXyxH/LM3oXplpUh5ejtN5/WmJhs1Ofb/QmH/tna31JmP7Y8ppvUnE8qIvZfvIRdrZp9MSNa+z/Qcaf7nfJp7k98mnWDcjjMi7ier2SnW1o51yhL7RS1vSpFclvefYsUB/BfDy1QxnBf83Rcu8aUSzjglIJGYIQuGfNgpCUlYScqWmAnJYxFMcBAEm6dZMeV242LgwL1iVOv+/5AQpop9WProO1/XW1WN+i/lMOqHyxzoQcQCApOwkhtcaXqD1199aj7+LP428Gr36YQuBQq7k2Nrb3DmnSablU8ONCjXd6Te1Ptu+vUileqVpN6KGjomiw8haHPhVZY7oMjawSBOr5b3x5c2ZkuP30OvzOuxYcNngOmkJ2ZSuoDKv/LP7Ib6BHmpnUIVcyaqJmtovY5e3Zu2MMFL0VMttO7wGzqVs9ebx6PRRAHuWXtVpLwzvftUIKxsLGvesyNm/dL8HnUtpNHlZaRot0rG1mhDnjh8GcOj3G9w9F0P/aRoHb3NzM+p2qICLpx37lpv23db1kyDMzc1Ijstk3YwwdfsHC4unmOp/FdEh1gQWD+uHLCuTyg2b0v1zw34ghSbyEiQ8hIC+fHDgA85Fq2zNFoLA5YinYOsKUyKKfl8RkUIiCAIHHh9g0nHjzuGNvRrza/tfTV636campEpVBdPCh4ZjJjHjUfIjuv/VHU87Tw73019duSgJP/KUU1vuabXldkBVKgUtzUD8szRS4jOpEOCOubkZCrmSrHQZ9s7WFDVXjz3l5KZ7+Q/Mh7HLW2tpVbqNC2L3Yl3fGu/qrjy9lajTnrOGMZ7eTsDC0pzt8y/q7e8zpR5PbyZoRUHl5b1vm6nfx7yOtAC1QsoR8m7VfB1XTXFQzc6QoZAL/D5ZtzBfDnmdd1PiM7mwN4IGXf1wdPtvZsXNj8Lev0WzTj7IsrOQZameGgLadCj6DZ6Hw6+tYNv7RPxQRS2YAJx4/LJezvj8HfJERF4nEomEjr4d8x13MUZ1Y3py8wUHV14nI0VqcOypyFNqwST33D9v/wlATEYMAasDCFgdwOnI02TKNU/zRWU+UiqUOoJJq8HVtCJj8posPMo7ULF2KXUosbmFWbEIJgABIeXzH6SHPpPraV2vmnRS61qfYAIUWjAB8K7mhlclba1K7bbeNO3tz4eLQyjj50yDLn588EMwY37W1ca1GlJN633UZzLJcYAdNtdwCnhTnVut7Syxc7LSMhOBKoS480cB9JhQRyeqyMnDltZDq4uCSTEgmnXyITZCI9U7lSraVMMA/KJRBfZ31dijN0RG45ij1LLRVZuKiJQExtYeq1O+IDcypYyjj49ya7Hq+t6FWEr7OtHvC+0QdX21fkYcGGFwXX1J4H5p+4tWFeuCkhSTodfRskbzsoVes6jRlwV1zM+tSIrN4M85/6CU6yrCcyJ0rO0syM5QmYzzRtMUhI/0CBLGGLU4hJOb71GtiZeOsAKo/XDGLm/N01sJOHnY4FxKN7QYoGxlF7Wj6gcLW6gFQgdXG0YvacnpbffJSMrmweU49ZzgdyrrXcsQlet74lrGHhCwdbDC3sVay3Qk8noQzTpGyEpL4/TmtVw5sAcohpo6Kc/hh2pIgX0O9kwvpalFce3RE9UP76yH6l31zxcRecMIgsCJZyeo5laNE5EnmBU2C4B5wfP44uQX2EodGXZxjt65ueuODNozSCsxYGEprE+UPpOBW1l73pnesMRV2s1twnh/QTA29qqbu6AU+DlXltJmff2p0awsVi+zlBoqOpiXTqMCcPSwYfP/6TqmwpvN4aGQKzm56S7Vm5bF08/wvSAni29R5XURKTxinpMiJiMlmWUjBxXvJpdUWWfr5UlNXys7VwhctS7FewYRkVdAIpEQ4h0CQL8q/ehXRRXuLlVImX14LkMuzTI4N8eBs/TIlCIRTACUgrJQ5RP0FeQbMKNhiYy++GBhC1ZOOIFvoIdaMAGVVuW9b5shy1Lg4qmreTA3N6NJr0qE7XigbvOp4caTmwkEv1MZW0crKtbRmKdyaNjND99ADw6uvEH793UzSb9OzC3MaDko/0hBUSh5+xGFEwM8u6n9BFalcXMDIwuJPJuHp7+jh56aOWuiYqBsXRh5FErgl6OISG5ObLpL/NNUeoyvQ0aqFHsXa6zMrWj/wLQondhfnaCJ6ue6pevydZOv6bGzh864a8Ou8TDpIQsvLiTEO4SlV5YSnxmvNSZoTRCXBl9iQugEjj87zsr2KwksFaiTmwXg0sHHXD8eSd8pulmQzcwkJVIwAVW9FkPaC3tnazBiBa7boQKPr78g6l4SA2c2emm+0E/ePQb9r3GhzisiUhhEs44B7p8/y87vNepo7xoB9P96btFtMNOZ6R5u7HTUzfJ4reoYqDMYrB2Lbj8RkWIgK03GbxNP6rQ36OLL+T0RWm3Lm4wHQcJ3kjU8DNMuFpdg+5zNteepI3RkShkWEgsC1wRiobBCKZETPlzbabMgtZ8m1p+I8xV/Ik+alpjrnekN1BlNRV4viVmJtNik8sWbGzyXrhVFs/bbjGjWKWLMLbTfmjL+upV6C4VcijCnFP/YWOsIJkNqDGFyg8lFs4+ISDFz/fgzjm+8q7cvr2Cyu8ZLp1mJwFKn6ThX8KPpY01NJ7dMLyzl1mqTjKWZylzR2bM7Pn+1AeBeQAxJMRk06OIHwI+tfmT8sfEAOFo5akX65GXNsS30vv5Zvq/po59blTgfk/8aOYIJwNSTU0Xh5D+KqDkxwONrV9g6Z7r6ulGvd2g+YEjhF5RL4eZfCNtHEuJTjkRzTWhi7VK1WdF+hV7Vs4jIm2Lt9DOkxGfh5e9M548CSYhKY/ficEpVcOT5/WST1zlU+Q8eeOhPCjY67Eeta3NLM3wDPChXxYWAluWN5q/o/FEAyV6ReDt6427rblSTkncffWRapOIyKo4BVQfgYPXfqSbb/M/mJGerfp97eu3h74d/k5ydzIeBH+Juq3HS33xnM5vvbGZd53XYWBR96GxydjLN/9Q1n/ep3IeZTWcW+X4irwcxz0kRkxIfq3Xt5FHKwEgT2dCfhL9GUdfXW0swAVjVYZUomIiUKCLvJJISr6qj8vx+Mr99fpIdCy4jlykNCiblq7nqdcQ0JJgAKEfd0LpWyJQ8uBTLiT/v5ptYa++ya9QuXVt9A50fMl/vON8X+lPMJ1vHaV2vbjCdHy/9SJONTdh2dxsAk45Pot/ufvqmFxvxmfGceHYCpZGs1EceHyFgdQAnn+ma1ArCsH3D1IIJQJcdXVgWvowNtzdoldGQK+XMPjubO4l3mHh8IoIgqHPOdNvRjXRZOjJF4cOTAX67/pve9m33thXIhAeozzb2yNi3uoTCfxnRrGOAg8sXa12bWxauJgZKJTwJI/VRKCG+uqXlzw86j6V5IdcWESkmXkSlF2h8m/eqU62xl/o6M1XK09sJVK7vifudL1l8eTHzgucx9shYrXnv13of6SxVZdrCkJUuU0esdPTtSEffjgiCgEQiYd3Ndfz98G+Cw95Xj/8z6BuS7GLU145Zbgy6/DVnKuzQWnfx5cX4OfuxP2I/AKcjT9OsXLNCnbEgZMgyaLVZO4/IpSGX1GauHD4N/RSAMUfGcGXIFczNtB94DBGdHk1pu9Jq89mlWP31Y3IIfRpKS++WXI7VCJjHnx0ncI1G4ItIiaDxBpWz7NZuWyljXwZn64LlZpIpZfx+/XejY+Iy4ihlp/8hMUueRWxGLF72Xur3BuDEsxPUXlubq0OvIlfKqbtOVVV+f5/9lHNQJWc7dS+eiqXsKeuiekDMmwFY5M0gmnUMsOAdbTtn10+nULVJcMEWUSpgQTVIjyVAT1ROjvOfiEhJIz+tRQ5N+/gT0LKcVgZVY+R2dgRNXpKDv93g3vkYQ9P4aGlLstLlnNl+nztno7X6hsxpgpWtBda2FlpJytISs7l3IYYz2+6r27JHXmFC3QkkZSepz2FrYauVbVYfjb0aY21uzfFnx9XnfpT8iAF/D2Bl+5UElCrYk31ukrOTkUgkNNuoX/gpa1+WA30PqK9HHBjB+WjtHCSbu26m/9/91dfH+h/Dw9ZDa8zKayv58ZLGvHWo7yHabW2X7/l2dN9Br1298h2Xm2+Dv6Vzxc4mj9/9YLdWJfbSdqXZ3n27lplnVYdVNCjTQGduhlROo411CnQ+gNRb84z2R8wrGWkcFEoFs8/OpknZJnTwLYYs5cVMYe/fonCiB4VczqJBPbXaek6eQaV6BShEdms3bBqsvswtnPg5+7Gxy0bsLQ2H8YmIvE5iH6ewZe4FAEbMb86qSdr1RfyCPGja2x8XTzuOb7iDpY05TXv7F2qvqLQojj87TluftlpPwoJSIP5Zmk5VWgdXa4bN1dy49SVMA3AtY8fAmY2Je5qKo5sNv32ubfJ4Z3pDPMprfEkUSgVyQY61uTXZimzqr9MNKTbEvt776LS9k/q6sMnf2m1tR3R6dP4DgVMDTmFnYad++s+PvGcyZhpZ1WGV0Yy8heHi4ItYmVsBqt/Z1ntbCfAI0FvRvNH6RmTIMwD4o+Mf1PNUpduPzYilzRaVQ/TiVou1ik4mZ8gImnUQUOBY/csCny/90ccoswyXAzjyeQiVSr1536PA1YEIqG7TpwacKrBW6k0jCidFyK1Toez96Xuttj7TZuEbZNqXAgAzVR+gpS7OLHfVfJjmt5hPR7/8a5KIiLxODGlKOo8JxC/QQ29fcSMIAmmJ2Ti4WuvkHHly4wW7fypYzan8MpsmZydzOfYynxz9pMBnvTT4kknmWYVSoWWCKagvRUHIq5k1tte1Ydc4H32eU5GnGB00mqWXl7L65upXPkOOgJR77+VtlzP+2HiyFaqw7sWtFjPu2DidOfrOndO36PBdFh2+C0gwt7uPXYWVOnvXdK/JjRc3dNpzyIruiizRcP6qOj4uTOpQFXd7a6qWeXNh5blfv6OVI2fePfPGzlIYRIfYIiR0je4HXWJWgLcq9pb6x9yCCUCbCm0KfS4RkeLAWBXWNyWYgCr7rKObjd5kaD413ek9qZ6eWfoZMqdJvmOcrZ1p6d2S8KEaoaeRl2na0rrr6pIlz9LblynPVDto1l5bm4DVASRmJerUEsrN1aFXCe0fanTPK0N06wDlJmhNkFoAiMuIMzju6yZfA9CgTAMm1JuArYUtExtM1Dv2UN9D6p8XtlzIuYHn+KrJVwbXznnduRl9eLT6XICWYFLRuaL659QsGW0WhGrN3fVgF/2Wn+G3yD44Vp+Krc8KrDyOqPtru7VU/7yyvfb3uK91MGn3NO+5lVsYAHfm6H9YvPwkiYG/nqPDohPsufrc4GssTk5Hnta6NhYu/29D1JzoIa+/CcCEjTsxM9HpTJjpzDZHe3xkct738tTqK6z6V0SkODDmW1KrRTlCBlZ9jacpOKb4xrQcVJWawaZVps0hQ5aBTCnD2dq5QNqNtZ3WUrt0bUCVSv+LE1+wL2KfSXNnNJ6Bq40rAR4BlLEvA8DO+zuZfnq63vHXhl1j3j/zWH9rvbrtSL8jajNIDu427thZ2vE09SkAweWCORmpMnmdfvc0Tlb6v1NTpak03agppLiw5ULaVmird2xkWiQdt6lu8s7WzloRQAXhypArJGXIqTfnsLrN3DYCO9/lJs1PvTWPW7M6MnfvTe7EpDG47QumnZqq7gNwrK4RUHK+j5MypEQmZeJkY8m8/bf1CiPF6YMS+jSUk89OMr3xdCQSCYIgkC5Lp8lG/UJ1jkDZt0rfYjtTUSGadYqQvMJJgQr+pUYTsF3XyczRypFt3bbh5eClZ5KIyOsnM03KqomGtSaDZzd5a6qxGhJShsxpgpPHq72G6/HXeXfPu5RzKMfOnjuZf34+m+5sMjg+xxl1853NzD47O9/1W5ZvyU9tfjI65nLsZYbuG6q+vjr0qpZG6YeLPyBTyJjScArrb61n3j+GnT0vDr7I79d/J7h8MDXca+R7PlPw/UJVHLVzQBn2XovWEgAKwuXB4fh/qSvMmbqePifX7/oGMnmrpnaTne9SzG1Vgtrytst1orD+uhzJp5t0tVLFJZwIgqAV/RT2bphBoSQvh/oeUguyJRVROClCCi2cKBUwy01vZE7eLxMRkTfNwytx7Fuu0eR9+GMIaYlZbJh5jvLVXOnxacEjIN4UmWlSts67oM7NAtBiQBUCWhp2eHwVpAoplmaWSCQSvZqVa8OumaxxMVWbmhPRMj9kPh19DfutZcgyaLRBvzmqbum6rO706r4kubkQkUDf5WFabbbeK7Fw0ERJyVOrk/lsECABzAEBJDIsnMKxLavKKZP5bCDyVP05aSQWiThU/tboObJiuiBLyD+i8vDEuvTa0159/Vm9zxheS1UHKikrieBNwQgKG9LuztSat3pEQ0Kq6IYyS+VK6s05xNhW/owOqZTv/lpnlmex6voqloUvMzpuUatFhMeF6w23HlN7DB8FfVSgfV8nonBShBRaOJnpzGZHB2Z7uGk1D681nM/q5Z86W0SkOLj7TzSHVt0EoNfndSlb2QWAw7/f5M45VaRIfs6ibwuPb7zAwcUaazsLHFyLPoupIR4kPaDnTk2EX0XnijxMfqh3rAQJAgJBpYJY0noJLjYuJu8jV8qxMDOenirvk3huzg08h52lbqK8wiKVK6kyXb/ZytLtBIoMX5RZug9r2ihRCS3aD29NKrqzbHBdXOys6LL4JDeiUrCruABz6zgU2aVZ3/N7tTZpfov5jP5FrrNGXs5/2ZZSjtZ0/6s7j5IfafXt7LFTp+DkO1UGsHKnykzn7WbLgQmNsbWw5VF8Oq0XHKd3nXJsvxypHj+zWw2eJmby26lH1K/gytaPmmKIgpgLcwTYb859w8bbG3X6d/bYibWFNSuurmBs7bGUtitt8trFjSicFCG5hZMh3y6mtG9FI6NfopDBbA+9WpOCJEkSEdFH5N1EXMvYY+dkVaB5Jzbd5dqxZzrtXT8O4u8lGsfPf4tw8iZJk6YZVMfn1Y6kSlNxtCq+CJCHyQ/ps7MPEokEmVKVuXVnz51aDqeFIUMqx0wiwcbSXG3KKQ78POw5+nmIlrZZEASeJWby+EUGDf3csLIwI12WjlQhxdXGlb+vRvHxBsPZiNtW92TlsPrqtQwJcHkRFNZkx7UDwQobr+1A/jlSQAGY6zUFZSuymXtuLtvubTO6gjShGdkx3QBVcfrxbSoTWN6ZJK7y9TnjD7uGzD2CICBXyl9r4k9ROCkibhw/wv6fF6qv89WaCAJEX4VfWvBpaQ+O2GueSnLH+Yu8OkqlwO2w55Sr4vrW+EIUBbn9Kcb83Eor0ZgxMlOlOvlK9FG6giP9puomtxIpOPqeht8GvwBQ3bgUSoF7sWkIAtQoq/1d++Pheyw8rCr0WK+CKxcfJ+a75sBGPpR2tGbR4XvqtoXvBDFhk0Ywvj27I9Vm7NeaV1j/joWH7vLjkXv8OKA2PWqX0xKgbs/uiE2uZIHPUp9p5aopKKl3vgKlthYqt29MZlRf5Mn1ufJVO77edYOdV6LYOLIxH55qqXe9jKfDsPPWmNxSb83FkCbo0OTq9N7dQ29fDlVdq3In8Q4LQhbQ3ldlxpoVNou9j/ayo/sOLf/HFSce8M3e2ywfXI+OtYr2syoKJ0VEbq1Jwx59CR74nu4guRQsXgodN3fB5iEM8vLkqo21eoiLtQsnB7xa3QsRbS7sfcS5XSpVbLO+/tRum5+6+O3GULKxVkOqUb2pl14fJkEQOLP9AX6B7uxYYPgpMjc9J9ShXFXXVz6viKouTu708yHlQ1jSZskbPJFp3I9Npe0PJ7TapnaqxqiXPhSmaEnCv2rP9J3Xmd6lOsfvxFHbx4UqntraoZzU8F/vvM7qsMd81yeQ/g28Sc2SETDzoHrcqzifyhRKLM1VqR9iU7I4eS+ebkFlsbLQTQdx68UtErMSGXV4lLqtf5X+PE9/ro5oMkZ2XGvkaTVRZpVDYpGMQ+W5Wv25NSwWjlexdL6MheMtrTFp9ychyB1BsMLc4SY2ZXaS/nC8juCTl0dzO5us/QFY3XE1w/YPA0Ca2JBFbebQKcBLyzT3fnM/ZnQtGifpHEThpIhYOmIAWelpAAyeuwjPinmyYN47BOtfhm91/h72TkQAAnOZc96WL6S3jbwRGSXNFCEIAinxmTh52OoXHJSCyVoPUOUfyUiR6u2r3MCT9u/X1Gk/vvEO149H6rSPWdYKiUSiV5tS0t7H/zJyhRL/L/fh52HPsYktDY4TBIGULDnOtqap5wVBIC1bjoO1hdZn83lyJk3mGg7HjpjXBUEQ8Ju61+j6Y1pWYnJH3cyvxsiWK7C20GgyWs4/RsSLDC3zy+ti692t/C/sf4DGBJdX0CwM6Q8/RZldBpty67F00nV8Tn80FmWWpubaqSmteBiXztTt15jWuTpRSZn8395bOvNyOD+jKa23aspBWJpZqs14xhAECWm35/Jobmet3+2mDxvTqKK7kZkFp7D3b7HwXx7qdOpG2NaNWNna6QomAOv7ctPKEielkrJ7JyJBWzAB+Km18bDA/zJymQJzCzMkEgnH1t3GzEyik0sjLTGb22HP8a9fGqVCwM3LnqSYDJ21lAolZuYlJ49gXi1HyMCq1Gqhyq/x28STZKXJ8PRzou8UzRfv9u8v8vx+MsPmNtVy4Lx+ItKgYAJw73wMLQdWxcpW8yecEJWuVzAZ+k1T9Q3J1tGKUYtD+GWcqkZMp9HFl6FUxHQypQpsrcyZuEVl7ngUn05KlgwnG/3CR9efTnEjKgUAMwl4Otmwf3wLnO1U4+PTsqk/5zB+Hvas+6ARzeZpBJBDE1rgZm/FoJXnuB1tPKmX7xd7GNfacJmC6//rgIN14W4juQUTgNBJrUjKkOJi9/pN4X0q96GcQzmquFZRt3nYenBt2DXuJ97Xqi2U9bwnNl5/mbSufcVFZMe31CuY5PVbcbO3oryrHeVd7Tj9heaBYWSLilSdvo9suW6V6gazz7B73BEGHlLltjk36Bz3E+9r1VnSh0Qi4FB1Bp9v0c7/U9SCyasgak7yELZtI2c2ryewTUfaffixdmfSE8780oBRZQx7Qv+XnF9lUgUX9jzi0oEnlKnoTJ/Juhk7pVlyfv30BL4vM41GXI1HYiahWR9/Tm3R2KGDWnvTrJ8/qS+yWDs9TGcdY1QIcKe0jyN3z8fQZUwgLp52rz1s+865aA7/flOnfcyyVjpCS5Nelaje1Iuwvx5w67Qm2VOOBkOWrWDF+ONac0b9FMIvn2i3AdTtWIEmPVWqd0O5PkTNSMnk2O1YHsans//6c85H6PpvVPF04OCEEOQKJaF34gjydqGUozVPEzII/k7X3AfQr155FILA9ku6QmpRcHxSSyq4//dqgsVmxLLgwgK6VuxKcPlgo+Ha0xpN45tz3+S7Zl7hxJgpK1uuYMZf1xnW1Jcui3X9yPLOjUqLYuPtjZSyLcX8C/MNrpvxZAQS8wzkKUFM7liNMS0LVy/LGKJZp4g4s2UDYVs3ENSuE20/0C7v/sOPFfjdxfD+weWC+bntz8V6vhwEmYz0f/7Brn59zKyt859QhChkSpZ/EqrTnmM6yI2p1W2Lg+HfNefWmSiy0+XU7+KLlU3hnvCkmXJOb79Pwy5+2Lvof6+L6nXW61SBGyejyErTqGZz3leFQsn5vx9xcd9jrTmNelSkfidfvWfoM7keZSq+XYXC/iuY4scxvJkvv5+OUF8/+KYzlaYZN7G8Co/mdkYpYHCPklKpt6TwPO057bepnE3tLe2ZFzyPlt4t6bqjK49THuuMN5dYoBDkZEa+gzxFk0do7fsNCa6sm0NFHy/SsrUy6AKse78RzSsbLjXx0aGPOBV1irS703GoMkfvmOLKxSUKJ0XE6c3rObttI0Htu9D2fU1iG+mZJdS794vBeb+0+4VGZRoZ1ZookpNJ/HMTNtWqYt+0KRJLlQpW/uIFsQt+wKlTRxTV6uPgYo25pba5Qp6YSMzsOTj36olt7drcbdBQ3ecxdiyYSfAYM+a1aAzuXYjh4ErDBbUGz26Mcyk7MlKkRuu2FJYBXzXkz1n/FGruoP81xsWzYHkefvkkFLlMpVLNK4AlxWSwZd4FpJnyQp0nPz76uRVmefxUzu95xD+7HxmYAeWquhLybhVk2QpKVyj+0HsR0xEEgUyZAjOJRCdC5U1y4NMWeLnYqM1It56n0OlHbYfQSR2qMrZV0T9Zv+0IgqD3ezdFmkKzjZrsszk5Zny/2E1OWTt7K3Ou/69Dgb+3M6UKFhy8w8pTmu+BEc38mN6lutb3RWqWjGk7rrM7PErddmFGM1pt1U1W171Sd/6v+f8V6BymIAonRcTpTWs5u30TtTt0pc2I0arGx2cICB2lNc7Owo5Ofp3Ydm8bzco2Y3k7w7Uf0s/9Q8LaNQjZUtJPav7gyy9fRuqhQyRvU8XOR/h04GHF7gBY21swYHpDzC3MuHTwCfb/G4y1NP96Ff6hx7AsU4bsTDnXjj0lsLW3QY2BoT8qQ+Q4dIbtuM+lA0+Mjh27vDW/fnocaZbC6LjA1uW5elQ3D4chRv/UEnNLM+KeprL5/86bPC83Qa29ad6/skljr5+I5PiGO1pt3cfXxru6KtFeXm1Fh5G18K9XWq8Wo9/U+myZe8HkcwaElKPFu/pr2xjT1BRGABMpfvouO8MFE8JvF71TW2/69LzsGx+Mp5MNWy485WFcOpsuPNUZ0zXQi79f1olZ/G4dxm3UjuCqVMqeb3oF6PU1OPMgnoG/nsPL2YaBDX34pI1pfzMi+XPkVgwn7sYxs3vNQj9QGnJU3vRhY/55lMCCQ3f1zouY14VBewZxNf6qTt8XDb9gUPVBhTqPIUThpIjYv2wRN0IPE9i2C3cvqm4MHcv3Y0ipAFKtX4AELg+5rM7SKH/xgoWWqTUAACAASURBVOx797Br1EjrQ3ZhXwTndj6ka4VwMlavMGnvoy2XGu13SH1Kw4v5Jf8B//DrWj4LuZ/2o6Z9SfL27eo+Cy8vKh8zbpJQKpQsGxuqt6/DyFpE3k3UccTM62sR1Mab8CPaX54NuvrRsKsfyXEZXDrwhIZd/Ti15R7+9Uqzf8V1ALqNC8KnhmEnrRz/jMDW5ZFmyrkdFm30teSl7fAaVG2kP65fn2CS+/Wd3xPB+b8f6bRLJBIehcexd5nGCa73xLp4+buwdvoZdYp1Bzdr3v2qEZbW5iRGZ7Dxf+e01srPV8RQqLHoY/L6yZYrsDQz09Fy5bD5wlOt+i76+HFAbcq72lLH25WqM/YhU2i+mrsHlWVXrqffMk42nJ2mW+H8h4N3uPIsmY9b+VPHxwVLczNiU7O49TyVkCqlSEiXUne2qrJw+NftTY72ESmZTNwSztaLpj/cfdc3kP71vQ0mDJwXPI8uFYvWdCcKJ0VE7jwnNq66WfiqDsmgyupdOHfrzouVK8m+o7l5+W3fhk2NGjo+GfUvfovMwp4X7jWocl83K6AAyC3sONncsONSDpXvbaZ03GWspSm4j/yAF7+u1BmTV8hxdLOhauMyBDVx4VFwM53x1W8bDlUDWDXpJJmpuuFpto6WjJivUg8+uBSrFijy4hvoQZcxqnj8zDQplw48oX5nX6xtjfuA5OREKAwR1+KxsDLHo7wDD6/EcS30GfFP0/SObTu8htqZ9f3vg1EqhUKZo3K0JrlJT8omI1VKKW9VvofcAkWPCXUonye/iEKhJOZRCp5+TpibGIm0Zd4FYiNUkRsDvmqIe1mHAp9dpPAkZ8gImqXK0XFnTketKJTwp0n0WHra0FTa1fDk0M0Ygit7sPZ9jYNlUoaU2rNUQsTIYD++7FKDR/HptPo+FICH33Qu9N9GplRBlkyBq72YIPLfQO7PnzFqlXPi70805pyrcVeJSouio19HQp+GEpcZR0//nliaFa3AKgonRUR+wkluXJLuUufKj0iAdNvSXKw7ieq17ZFZOnL7SoreOWYKKcqXWWN9nhzmiY9uCfLqt9dyq9oQo3vXubKIJmf+4vnUqVj5+5N96zbJe/dyv1JvnnrrPlHlUPXuRspFaW68OYJMQKvyNO/rrxOaG3U/iR3fX9K7Vl7zwYuoNL2+IKOXtjT5RlucJESls3HWufwH6sGnpjtPbrww2F8QbUVGipTUhCw8fYvus5yVLsPKxrxEhVaXVNKy5Ww+/5R3G/pga2VaZJ1UruTm8xQCyjljbiYhOVNG0P/03xAGNvLhm16qEO38nF4fze1cILW+QimgFAR1kjEREVCVFZj+13WDUVp/ftiYxm8oTFgUToqAxOh0Vo1/R32dn3CSg1V2ElJr04t35UfrUO0oIUPmntzmmtTjx1mz0bh/R971H/p2IcK3s07/iPnNUcgF7JytWDZGf8hil7GB+AZoe4cbMv+UJDNDRoqULfPOk5aQXaB5OY6py8YcQ6nU/pPp82ktPKuWEqtOl3DuxaTSbqF2FlRTo09yCxkXp7fViZbIy505HZEg0SmKd2lGO9zsrTh9P556FVy10qmLiBQFK08+xNXOivi0bKp5OWFlbkaTSm8uf4konBSSxOh0Nsw8h2+gBxFX48lK/OFljyU2rp8U6V6mULNFWZo0MCdh3XrSjh2j1KefklW7Ndu+u2hwzpifW7HtuwvERGgSKnnEhxNwfQXH9Ag2laIP4zfuPQ5vNt1WCTByUQvin6apTA56UkGDrh/EiO+bY+tQMtXHD6/EsW+54XL1ZUpDdKzq5xwBSxAErh57xqnN96jezItAi2vE/e9rAPx27cSmShVDy4kUEyfvxTFl61WikrOoVMqew5+F6AiKh27GMHKNrjNy+FftcbazJEumYO7eW3SvXY56FXRT+RdFkbt944Op7iVGT4n8t3htwsmJEyeYP38+Fy9e5Pnz5+zYsYOePXsanXP8+HE+++wzbty4QdmyZZk8eTKjR482ec/iFE7yRj3kCCe3KiRjJ0vAUzmJCrEy6lxZhNTKiQTXatzOx+SSQ/8J1cHWgc3fGI4q8fRzovv42ihkSjJSpLiVtTf4BJ6RIuXBpVhO/KnfCzs3nV4sJ/ua6sab4FoVpZklVwM+ymeWcUZOr8LjIUNxGzoU9+HvGRyXlpjNuZ0PaP5OlXz9SooTZXo6EisrMDNDYq77hCooBX5+qRn6aGkI9zt04ZncC9vMOCp/PoK4r75Uj/WaMxtBLsd1wAAEQSB52zbsGjbkQfsOWmvm578D8GzceFIPHqTy6VNYuJecjIxFzc4rkYz/8wq7Pm5GYPlX0yzKFUpuPU8lOVNGM393rb8RfYJDjho7OVOGk41FvunX85JbozJ2wyX2XH1uZLSKvKnADa0nIvJf4rUJJ/v27eP06dPUrVuXPn365CucPHr0iFq1ajFy5EhGjRrF6dOnGTNmDBs3bqRPnz4m7flahZOk5SBk0Ph+JG7pWTrjq926ycV9ESRGZ+Bdw40Xz9J4GB6PnaMV1Zt6cWzdbVw9banbyZdqjb105hfHmfOSU7lW9vw591tpTCqGzEN1L/+AS/IDrtd4n9jSdXX661/8DtfS1igeaISiinv+xrpSpUK+AsMoMzKQ2KoqDsujosgMD8e2Xj0sPT0LtM6L31YRO1/XwdjM3p4qF86rb25Px4wl7ajpCdTKL13Cs7EfG+yvdPgQVuXLI4+Px9zNDYmZRsMkKBRIzM25Va26ui0n9PvfSG6hoWugF5VKOTChXcE1SwsO3uGno/e12la9V5/W1TypM+sgiRn51xLJTetqpTl6OzbfceFft+f7A3dYe1Y3mRbAoEY+jGpRicVH7zGudWV83O1YFvqAb/ff1hkrCici/1XeiFlHIpHkK5xMmTKFXbt2ceuW5oly9OjRhIeHExZmWpry4hJOLu24SNgB7dwhqSnfY6kwo8m9Z7hmaPslmLu6UiXsTJHtX1geX3/B30vC9fYZ8u9QSqXEPsvUMQ/5BrrTyOoicd+rbuRyc2se+XbF+9lRnpULoVzUSWyzdB1By3z9Fa7vvmvwjEqpFImZGRILbc1J9P99g021aijT04n55ht8t2zBNqAWANmPHvGwk64PDGg0EpGffUbK3n1Y+fpSaf8+vWMBLQHgdeI59Qukjx+TuGGjus11yBAS1641Oq/ajet6tTtvKyfuxjF0la5zdDkXW626IbkRBIHNF54yZZtK4+fnYc/2j5pS52Xoa1GQIyRsOPeEaTsMm/T00btOOf6vVwDVv9pP/QqubP2oqd5xeTU5F6e3xd3h9WZxFhEpKZRY4aRFixbUqVOHH3/8Ud22Y8cO+vfvT0ZGBpaWumFL2dnZZGdrBIOUlBS8vb2LXDhZPPwnlIpYzC39sZPbEHh9BSf8VU/tuYUTm8BAKqxZjZmNjbHlXiuZqVKsbCy0QpZHL2lp0Bckh7w+IWOXtybr5k0e9TZNi5Ub14EDKfPVDJ12RVo6d+u/LG5nbk7VC+cxs7U1KDD4nzhO4voNvPjFcAbequFXQKHgTl1N/R5DJhR5YiL3mui/cZRYLCyofr1gN8uSSu4wWH0YihyYuv0aG/8xntzPEJ+1q0Ljiu7M3HWDm8/1R8qBtgYjtxAxq0dN+tf3Zvjv5wl7qD8q68E3nTE3k6BQCphJMGh+3X89mtHrLuJsa8mVr9qJjtIi/2kKK5wUezxadHQ0nnlU8p6ensjlcuLj4/XOmTt3Ls7Ozup/3t7eese9KgpZBIrsyygVcTQ9NxOH9Cgs5KqIFwuJGdVuXKf67Vv4bd5UogQTUFWXNbc0Y9TiELp9EsSoxSH5Ciag+kIdPLsJ1Zp6MWyuKueJdbWClTrPIXHDBm7XqYsiVeOIK33yRCOYgEqgqFOX9H8Mp5u/3yLEqGACcCeotpZgAirtTF7uNm6iJZh4zZkNgOtQ0/yEvH9dgZm9bmGz0pMnG5zjt/Mv7Bo3Nml9g8iLJ/3960auUBoVTAAGrDjL/uvR+H6xhx9eZrG8F5Oar2Dy4JvO/Digtt6+cW0q09DPjb3jg1k2SNc0CehobD5qqTFLDm3ii42lORs/bMykDrpZeT9sURHzl3lFzM0kRgWOjrXKEDGvC+FftxcFExGRQlLsmpMqVaowfPhwpk6dqm47ffo0zZs35/nz55TRY29/XZqTNZ07EedoTuXoNCrHxABwoJYfCnMzetZtTqUpXxTZXm8DiuRkMDNT1+1xHTgQl759tLQqPn/8zpP3huvMrXY1HKVUyt36DV75HA6tWyOLfo7Pr79yr1lzo2Pdhg3DtnYQFp6eKBISePaxdoRVXu1K0rZtyGJiyLp5k7TDR6i4dw/y6GiejHhfy8SUdesWCevW4TVrltrcklvzU+3WTRI3bEBiaYlr//46/j0A5RYtJPLTCXrP7Tp0CHb16xM5bry6zalzZ8r9sCCfd6dkM2XrVa006n9/0pxa5ZxRKgUqGigmd/izFrT94YTevhzGt6ms9lcRBIGfQx8w/4AqAeLPg+rSOUDbv0uuUKoFg+iULMq52OqsmZIl47v9t+keVI6Gfm5afcmZMrJlCko7layHEhGRt41/lVknL8XlczJj4lBcniZwsUoC1/xTQRAYvrcCgkRClfHf0a1pDb3zsmQKBAEikzKpVMpwdM3biiCXg0Si5QMhffwYc3cPzB3seTzsPTLO6SYzk1hZIejRZuSl9OTJpJ86hfTJE2TPtMOZ/U8cx7K0Jstq7Pff82Llb4V6HdVu3tBySH1VBIWC9DNh2Ddvpvd3/uT9D0g/fRrfLVuwqVlDvXfKwYNIrKywrlSJB+3a49CmDd5Ll6jn5TV3lf/5Zxxbt9J7hpS9e4n87HN8N2/CNjCwyF6bMQRB4OONl7UiViLmdeHXEw8Jf5bET+/WQSKR6NT6GNOyEpM7arRy+64956P1+hP65ebrbjVo5u9B+1w5SfQ5lP4cep+opExm96j1r/sbFBH5t1DY+3exx3k2adKE3bt3a7UdPHiQ+vXrmySYFCeNfAK48/Q4Uqk3cBN5WhUEiermOnnHDSyd3OhYS6PZUSgFvaXE787phJUJJpW3hbxOrABWFSqof66w+g8Ann3yCamHNMmo8gom1W7e4HaNmlptlfarnFndRwxHEATiFi7CurI/Tp066d231OefI8hkJKxeA6h8T2RRUQYdZ3PvrU8wEQQBhVLAohAZNiXm5jgEG9bk+PymKiWQli1n2l83GNTIhxtRyfhXrU+9Cqonc31+MnaNGmkJe8/GjKHK2TDMXbTDb5+M/FBdODKi/zsmhS2/KhlSOe0XnuBZYqZWe25/DU8nG2Z0rUHgTE3G1CUD69A1sKzWnE4BXqwYUo8P1xrO2bP+g0Y081cl9+seVJZ9159za1ZHvWPHtBQr5IqI/FspsOYkLS2N+/dVYX116tThhx9+oFWrVri5ueHj48PUqVOJjIxkzRrVzSQnlHjUqFGMHDmSsLAwRo8eXSJCiff+9D23ToVSrccgAjp049uJU6iYoQobXOU9hHQLVY2SiHld2H7pGZ9t1h8hA+BoY0G3oLLqtNX/FQw5uVY+dRILDw+yHz4k7cQJ7Js0waaq/gq7hSF+xa/E/fCD/k5LS6pf0y2ydvVZEt2XqOqcjG1ViUkdCudrk8O1Z8l0W6IpBRBQzpmA8s5sOKfrO7F3XDA1yhr+7Op7Hz2nTSP+l19QvHgBZmagVGr1l1++DMeWLQv/Aoyw9eIzJm4x/HnPjxznUX2sO/uY6X/pr8MkhtyKiPy7eG1mndDQUFq10lU5Dxs2jD/++IP33nuPiIgIQkND1X3Hjx9nwoQJ6iRsU6ZMKRFJ2HKEk5DBI0hLfMHFPTvVfcsqfID8ZQGkB9901qsx0cfIYD8mdajG7egUapZ1NvgF/W8i89o1Ivr1V19XuXABcwddp9Ki5MitGMqGh1G5QwhJW7cSt0AlqPgfO4qll/78MnlDPAc18mFCuyp4OFhz8WUpe33ZQfMiCAKn779g8G8Fq9Mzt3cA7zb00WlPy5bD+XM8/eCDAq0H4D5qFEJWJh6ffIK5g0qYzr5/n+jZcyg98XNsAwJI/HMT0TNn4tKvL16zZ+td50FcGm0WHOeT1v583r6q3sRmjtYWZMoUyJXGvzJqeDmxd3yw0TGgKuy46cJTqng6UMfbtdCF7EREREouYvr6QrB3yQJunTxGxboNeHhJO4vrEt9RCBL9qv+17zck4kUGTjYWjP/zikl77f64OTXKOv0nhJXiILej5ZKBdfh4w2VAFYGR4+woCIJe34OCph7/a2wz/Dzs+ftqFMH+pRj021meJqjMGnfmdKTq9P2Ffh15C73lnO3DFhWZ1rk693v3QXbzZqHWrrBhA0jg8bsD1W0SGxuELE0ywWq3bqr8Q5RKEn7/A/tmTblp50nvn43n7+lbrzzf9wsCoPWCUB7GpRsc+yoVc0VERP5diMJJIcgRTvTx2cZd+E3TTfI1u2cthjTW+F9ce5bMj0fucfhWjEl7ftm5OiOa+xnNkyCi4kFcGnGp2QxYcdbouKDyzuz8WNcXJEMqp8ZXB4rreGrqVXBVa15ycLKxICVLNzx42aC6uNhZIVcqGfKbdnj1ltFN6Lc8DEuFjF27p+rMBagcdgZBKuN+SEihzmpdozoVt28ncvJkUnapfME69fze4PjZPWvRvoYnnkaiVm49T6HTjypfmC87V2dki4qFOpuIiMi/D1E4KQT7lizgphHh5GliFi3ma/cbsomfeRDP/7d353FR1fv/wF+zMcM6CsjAgCIom6JGoKGi5IqJVtffdct9Kbk3ErRM07qZqdDXMtPU0lt601wqLJdMRTOXQCkFRcF9gVhElFWYGWA+vz+II8PMwDBsA7yfj8c8Hsw5n3PO57xHmTef81le2Vq/Zv4qqStGGbx0e3uRVVCK/lGGTyt/8q3n4WZviQv38/D/NtfeCrBjTj+txMBYY/vIsWGyHwDgSkblbMPFynKNScZqjmIxhKiiDN3zM+BQmodTzs8AfyeyEcM88PqQ7rjt27OOMxjuvf5z8KdMd98h6gNCCGkISk6M8MvGtUg5rfsLcOGeg+DxeEh7VILNp25h/jAPOEm150qoLu+JClYSIUQCPq5kFGDMhrO1lq+urX4JfHgoBV+dvQsAWBTihTVHr2PthD54WKTE814O8HK01ij/3OrjsLMU1zrLpy5ScxEKSuteY2XmgK5Y/mJP7Dx3H9kFCnSxtcDbMdqdZ+vyT38XLB3tA1tLw1ZcboxVbav0cpbiS6dcFLxb2brievo07g8erFFmp9cITL3+dDI05uQMXlaG3nO+FfRvrPeugOq/X+J4Z3984j8ZsQsGw0NmrfcYQgipCyUnRjAkOWkIxhiU5WpIRAJk5pfi7R8u4+wt3bPizh/aHQtHNt5olpZ040GRxhwVtbkbNRq3Hz7BikMpOH3jYa1l70WH4uaDImyLu4d/+DlDZi3RatnSZVRPRyx/sSccpbofTRxOzkL64xLMCXKD7/KjUJQ9HRVTPWl8Y3ciDl7KxH/G9MDsIDeD7q/KgUuZmL87UWv7iTeDYW8lRp8Png7DTVg6DBWMQWougpoBO+Lvay8mxxgCcq4hw7ITsqzsIS5X4X/HVqHIzAKvDl8MADAvU2BgZjJOuzwDlUCEX356y+D6vr/ga3w3r3+97pEQQmqi5MQIRzZ9iqunTmhtHzo7DH4hYxrtOtUlpuVhx7n7eDe0BzLySjWGogKtvwWloLRM44u2Lttn9cXMbX/o3Z/0nxEoKC2Dq53u0T//2X8F38TrXjUWeNrRtD4y8ktxKT0fL/g6Nlq/oOozpC4K8UL/bnbo4WQDiciwx3kpmYUYvf6Mwde7s3q09oysjGHLif+D05NH2Os5DLu8huPnA4v1nsPzzz+4EUCEEGIMSk6McGTTOlw9dVxr+5t7DzXaNepSc9bM8CHd8ZaOtT1MHWMMl/4qwMsbf2+0cy4Y7omI4R51Xrd6fw59i8qZAlW5GhVqZnT/IlW5Gp7v6l+JGQA+ndgHwZ4O3OMmn/eOoLSsgtsf1N0eVzILkF9ShnnB7vBT5sB1if5h/c0x0Zvy9m0U/3YK5TkP8Ph/38Aj7ncIbW3rPpAQYvJMdoZYUrsXejlhRA8ZYlMqR/t8fvIWXvaTo7uD6T/rzylU4NDlLEzv74rpXycg7rbu1VzPvTMMRYoy5JWUwUosxMnrOdy6KNUFuHZE1LhesLMSI6ugFD3l0jrrwOPxkLpiFHacu4dxz7rA3oSXpm/oLMJmQj6OLRis9cjMRiJE0n9G6hy+m/rhKOxJSMOSfZUrHr850hN+XZ7O5cLUXrhWyxJSZTk5GssJ6FKelwe+RAKeWMzNCNx561aIu7mDqdUwc3HRe2zplau4989/amy7OWAg+JaWcF7/GawGDqz12oSQtql9t5xsXoerv7Vsy0mViD2J2J+Uyb2vOR+GqaltITegMiHJLVbC11l3glHz8c+0QFd8+LJvo9ezLer1/lEUKcsxsLsddsx+zqA5RdRqhiJFOaQWupeMUN27B8W1axB26oT7U6Zq7Ou6dw8kvXtDXVzMLQop/+RjSENDcW/iJJReqn0mWYvAQLhu38a9L8/NhSo9XWM+Fn0cP1yBjuPH11mOEGKa6LGOEY5+8RmunNRc3v2fy1bCtbfuZdmbUlmFGh7LtJvsX3pGjs8m+TV7ffSpUDNs+PUm1h2/WWs5Q/vOZBcocP1BEQZ72Jt0Mtbe6FuWoCG8L1/CzcHBqMjPr9dxuh4tKW/eBAQCiN1pThVCTBk91mkkLZGYAIBIwMeC4Z749PgNje37kzLxgq+TxgKELUWtZ+HD6urboddRKtE7ioa0HPdfDte5uGJ9XevdR+8+53XrYDMqBACQv+9HZC1dyu2rSpR45uZgpZoLEDosegt2c+bUet3yhw9RfOYsKvIew3bmTI3Vtgkhpqmdt5ysx5WTmiNLWuKRTnWPn6jw7IexWtuvrxwFsbBpfqlWqBkupuXh2S4d9U6v/6hYCf+V2o/Adszph3k7LqBEVYFFIV54fQitFNtWsIoKXOtp2KM27+TLAJ+v9cWfuXgJCvbv13NUJcugIHT571bNazOGaz49DLq2/OOPIR3zNCkuPnMG6a++ZtCx1TlvWA+bESPqfRwhRD96rGOEY1+uR/KvppWcVBn92RmNicgGe3aCh4MVZg3sCsaAzrYWDb7GzQdFyCspw4Qv47ltVZOUValt8rAOFiIk/WckGGMoVzOIBA3r8ElMT+mlS7g/fQbcD+zH7ZBRAIDup09Bef0GHm/bBp5EAsdlSyFydtZ7Dl2PiKyGDIFFv36wmzVT73G3Q8dAdfu2QfW0HjEC9uHhyPn4Yzw5Y/iQ65rM/f0hj46CurQUEk9Po89DCKlEyYkRTDk5AYDyCjW66+iHUl1vFykWhXjBwkxQr5VdX9l6Tu/oGqByJebfb+Vi+tfa07xP7++KFS9R51VimCfnziFt5iwAgE1oKJw/0b+WT3Xlubm4GaR7dWP3QwdxZ8zYRqujLh0mT4LT++836TWaQ0XxE41VwplKhaITJ5D17ntwXrcOgg4dkB/zA2RLloAv0XzEWn0xTbVSCb7YdEfDEdNEyYkRjm3ZgOQT1RaG4/Hw5p6DjXb+xrD0x2TsOp9mcPnRvRwR4GqLFYdS0MdFivWT/XDochYupedj1T96QcjnIemvfMyqZeKz2qyf7IcX+8iNOpaQxlR6+TLuTZiod7/V8GEoPl45yaJ5gD8sAgKguHwZT+IqWwqlL72Igv0HIJTJUP5A98KdItcu6Hb4sMn3Uyk8fBhmbm6Q+PhUtmRmZiL3yy3I/+47jXJdtm/jEkV9rENCIFu2FOlzX4Xyxg2t/c0x9w1pOyg5MULsls9x+cQR7v2rG7+GjX3tczq0hMZcl0Wf/04PwPsHriIjv1Tn/l1zn0P/bnY0ooaYnPTwcC4JAQDpuHGQr14FoPIv/7KMDIicnWv9t3tn3DgoU/R/6ZryF/KT8wlImzEDACDo2BEVeXl1HNEw9uHh6BT+epNeg7QdNFqngSa+H22SiQlQOQKmSFGGsgqG5IwCzNDxqKUhxvk5Y3gPGYb3kOGvvBIEfaS5Xs2tVS9ASP1JiIly2bABWUveQXHc7+j2yxGNRxg8Hq/WSeCquH3/PRQpqeCZiQC1Gnf/MU5jf81HIzX3MZWyxWa1zdu162ld6pGY6Br9ZIjczz+n5IQ0uXadnChKnnA/l6mULViTullLKifPCvbsxA3XTX9cAiepBAI+D6t+TsV//179FwDEQj6U5Wqd56py7cNRWmu7uHS0wIcv9cR7+68CABLfG0GJCTFpPB4P8o+iG3YOgQDmvZ72o/K5loon8fFImzUbAHAjIABdtm+H4loqJF5eyFy2DG4xMXgSF4fMNysXVOwUGQn7sHkNqkd9qJVKKK6moOjoUb1leGIxHJcvR8FPP6Hk/Hluu9fFC+BbWEB55y7ydu9Gp/lvcBPs1eQRH4fiX39F1rJ3uW2p3j6we3Uu7MPCwLfUnbSpFQrwzMzA4zfs9wfXoVooBMrLAaEQXffsgbh7N60+MqTtaNePdQ6ujcKN85Vrwbz41jJ49G1bq7CqytX4OTkTQd07wUzAx7GUbAR7dYKDNf2HJsQQ9Z2MzvtKMnjC5vmbr7a6uWz8HFaDB4MnejojcHleHvhiMfgW+kf6lefl4Wb/Adx7113fwuLZZ+u8ptelJI3OstUTO0MfiT3ath3K69fhtPJD8IRClD9+jFtDhoIp9f/h2CkyAh3Gj8eDjz6CxMsLtrNmNTgZIo2L+pwYYePcV6AoqhyuO3bBEngGBjXauQkhrV9FcTFuBPSt1zFu+/dD4tX4w5DVSiWKjhyBzejRUCsUWi0dngnnIWiE349MrUbJH3/CvJevViJTfPZ3pM+d0zLJTgAAE3lJREFUq/M479QUsNJSXH/WX3tf8mWNRAmobFnJiFyA8txcSF9+GQ9WrgQAiH184Lj0HdyfNt2o+ktfegkF+/fD81w8BB06GHUO0ngoOTHCJxPHcD+HfbkDlh061lKaENIePUlIQNr0GfU+ziIwEPZh81B07BhsZ8yAmaurQcepFQrucUXavHl4cuo0OkyaiPw9e2s9rrk67aqVShT8+BOyly/X2G49YgSKYrUnkKzOM+E8wOej5Px5/PV6eBPWEpD07Am3mB+guncPQgeHWluM9FErFLj+TOXyId4pV6lVxgiUnBihenJiSvObEEJMS1lmJsoyMmDRt7IVperXZtUIoIKDh5C5aFGt55CvWQO1ohSFPx+G7J0lkHh5VZ5LpcKD6GiIOndBzkcfAQAs+vYFTyIxaEI557WfwGro0Bbpf1F45CgyIiOb9BodJkyA04oPuPesogLKW7dx96WX6jzW/t//Qu6mzQCALtu+hkVgIDfzsEd8HIQdtf8gVd69i7K/MlAUG6sxFFvk2gXda+nfQ3Sj5MQIlJwQQhoLU6txrUfPugv+Tf5/H8Fm7FiDp+nXxfOPBAisrY0+vjHUXAupilAmg/Ona3H/lSlGn1tfAqFLwf79yFy8pF7n97mWitwvvsDDdZ8ZVL6qI3F1ZZmZEEilejsGt3eUnBiBkhNCSGOrz7pA9dXlm/9xj5jcDx6A2MOjSa5TX9U7yjqtXo0O4/6hsb/sQQ5uBQdrHeedfBmKGzeQuXgxnNesgcTHB2WZmVA/eQKRqyv4ZmYNrk9jc1r5IWzGjMGTs2ehuH4duRs+B9C8naFbE0pOjEDJCSGkqSnv3MWd0cat8NzxlVdgN2c2Hm7aBPtXX4VZ166NW7lm9uTcOUCthuWAAXUXboDS5GTcGz+hwefpMH48bF4YhbTZta98XZ3nn39AYGXV4Gu3FTQJGyGEmCCxu5tWZ9WiEyc0OoTK16yB9bChlfOCCIVgKhV41VoN5KtWNVt9m5JlYGCzXMe8Vy94X0nG7dBQ8EQidDt0CKr0dAg7dQJfIsG1Z/zAFAquvNWwYRB36waHhQt0ns9m9GgUHj5s0LVvBPQ16RmFWwtqOfkbtZwQQpoTY4ybw4MmE2tejDHcf2UKShMT4XXhT4P6i9S3T5Fs6VLYTp/WkGq2CdRyQgghrQiPxwOPkpIWwePx0HX3rroLVj+Gz4dnwnmUJifDcsAA8Hi8yqULSkugvHsXiuQryFmzhiv/YPVqdJgwnhJPI7XrQdu9ho4EAPiHvtzCNSGEEGLqBDY2sBo4kBtCLrCyhLBTJ1j26we7ObO1yl9/xg85n6xt7mq2Ce36sY66ogIP7tyCzL07+Ca+JDohhBDTVpaRgVvDhuvcJ1+zBtKxY3Tua8uM/f5u1y0nfIEATh5elJgQQghpMJGzM3yupaLbce2ZcjMXLUJxtUn11CoV/oqIxOOd3zZnFVuNdt1yQgghhDSVmvO78MRiWA0ZgqIjR7TKmsKEek2BWk4IIYQQEyKSOcAr8SL3nv29eKMuN/r2w5PzCVClpzdX9UwaJSeEEEJIE+Gbm8Nt/36DyqbNmIHbI0Yi1dsHreChRpOi5IQQQghpQhIvT9jOmsW9d3jrTXjEx8HnWipc9QxpzliwEEBl35T8H36AWqVqlrqaCupzQgghhLSgkouJuP/KK1rbLfr2Rckff2hscz/8M8xcXaG4ehWq+/dh2b8/bgYN0ijT7cgvENjaQlDt+zLnk7V4tHUr+NbW6PbLYRQe/gUWz/XjVsduKrS2DiGEENLK5az9FI+2bGmUc9nNnQOHt95C9qrVyNuxQ2cZ2xnTYffaaxDa2TXKNWui5IQQQghpA5pyVWV9JD16oOv334HXyFNr0GgdQgghpA3wiPtd473X5UvwuZYK+3//W+8xXkmJ6P7rCYg9PY26piIlBdkrVxp1bFOglhNCCCGkFam5arUuNVtfnD9dC+tRowCAm37/8c5v8aBaQuJ+8ADEHh6NWld6rEMIIYQQAED5o0dImzkLkl69IFu6FAKr2ldeZmVl4IlEjV4PWpWYEEIIIQAAoZ0d3A8eMLh8UyQmDUF9TgghhBBiUig5IYQQQohJoeSEEEIIISbFqORk06ZNcHNzg0Qigb+/P85UWwa6pu3bt4PH42m9FAqF0ZUmhBBCSNtV7+Rk7969iIyMxLJly5CYmIhBgwbhhRdeQFpamt5jbGxskJWVpfGSSCQNqjghhBBC2qZ6Jydr167FnDlzMHfuXPj4+GDdunXo3LkzNm/erPcYHo8HR0dHjRchhBBCiC71Sk5UKhUuXLiAkSNHamwfOXIk4uLi9B5XXFwMV1dXuLi4YMyYMUhMTKz1OkqlEoWFhRovQgghhLQP9UpOcnNzUVFRAZlMprFdJpMhOztb5zHe3t7Yvn07Dhw4gN27d0MikWDgwIG4efOm3utERUVBKpVyr86dO9enmoQQQghpxYzqEFs19W0VxpjWtiqBgYGYOnUq+vTpg0GDBuG7776Dp6cnNmzYoPf877zzDgoKCrhXenq6MdUkhBBCSCtUrxli7e3tIRAItFpJcnJytFpT9OHz+ejbt2+tLSdisRhisbg+VSOEEEJIG1GvlhMzMzP4+/sjNjZWY3tsbCwGDBhg0DkYY0hKSoKTk1N9Lk0IIYSQdqLea+ssXLgQ06ZNQ0BAAPr3748tW7YgLS0NYWFhAIDp06fD2dkZUVFRAIAPPvgAgYGB8PDwQGFhIdavX4+kpCRs3Lixce+EEEIIIW1CvZOTiRMn4tGjR1ixYgWysrLg6+uLw4cPw9XVFQCQlpYGPv9pg0x+fj5ee+01ZGdnQyqVws/PD6dPn0a/fv0a7y4IIYQQ0mbwGGOspStRl4KCAnTo0AHp6en1WnKZEEIIIS2nsLAQnTt3Rn5+PqRSqcHH1bvlpCUUFRUBAA0pJoQQQlqhoqKieiUnraLlRK1WIzMzE9bW1nqHLBujKqOjFpmmR7FuHhTn5kFxbh4U5+bRlHFmjKGoqAhyuVyjy0ddWkXLCZ/Ph4uLS5Od38bGhv7hNxOKdfOgODcPinPzoDg3j6aKc31aTKoYNQkbIYQQQkhToeSEEEIIISZFsHz58uUtXYmWJBAI8Pzzz0MobBVPuFo1inXzoDg3D4pz86A4Nw9Ti3Or6BBLCCGEkPaDHusQQgghxKRQckIIIYQQk0LJCSGEEEJMCiUnhBBCCDEp7To52bRpE9zc3CCRSODv748zZ860dJVMRlRUFPr27Qtra2s4ODjg5ZdfxvXr1zXKKJVKvPHGG7C3t4elpSVefPFF/PXXXxpl0tLSMHbsWFhaWsLe3h7z58+HSqXSKHPq1Cn4+/tDIpHA3d0dX3zxhVZ92stnFRUVBR6Ph8jISG4bxblxZGRkYOrUqbCzs4OFhQWeeeYZXLhwgdvPGMPy5cshl8thbm6O559/HlevXtU4R15eHqZNmwapVAqpVIpp06YhPz9fo0xycjKCg4Nhbm4OZ2dnrFixAjXHHcTExKBHjx4Qi8Xo0aMHfvzxx6a78WZUXl6Od999F25ubjA3N4e7uztWrFgBtVrNlaE4G+f06dMYO3Ys5HI5eDwefvrpJ439phRXQ+pSJ9ZO7dmzh4lEIrZ161aWkpLCIiIimKWlJbt//35LV80khISEsG3btrErV66wpKQkFhoayrp06cKKi4u5MmFhYczZ2ZnFxsayixcvsiFDhrA+ffqw8vJyxhhj5eXlzNfXlw0ZMoRdvHiRxcbGMrlczsLDw7lz3Llzh1lYWLCIiAiWkpLCtm7dykQiEfvhhx+4Mu3ls0pISGBdu3ZlvXv3ZhEREdx2inPDPX78mLm6urKZM2ey8+fPs7t377Ljx4+zW7ducWWio6OZtbU1i4mJYcnJyWzixInMycmJFRYWcmVGjRrFfH19WVxcHIuLi2O+vr5szJgx3P6CggImk8nYpEmTWHJyMouJiWHW1tbs448/5srExcUxgUDAVq9ezVJTU9nq1auZUChk586da55gNKGVK1cyOzs7dujQIXb37l32/fffMysrK7Zu3TquDMXZOIcPH2bLli1jMTExDAD78ccfNfabUlwNqUtd2m1y0q9fPxYWFqaxzdvbmy1ZsqSFamTacnJyGAB26tQpxhhj+fn5TCQSsT179nBlMjIyGJ/PZ0eOHGGMVf5n4vP5LCMjgyuze/duJhaLWUFBAWOMsbfffpt5e3trXGvevHksMDCQe98ePquioiLm4eHBYmNjWXBwMJecUJwbx+LFi1lQUJDe/Wq1mjk6OrLo6Ghum0KhYFKplH3xxReMMcZSUlIYAI1fwvHx8QwAu3btGmOMsU2bNjGpVMoUCgVXJioqisnlcqZWqxljjE2YMIGNGjVK4/ohISFs0qRJDb/RFhYaGspmz56tsW3cuHFs6tSpjDGKc2OpmZyYUlwNqYsh2uVjHZVKhQsXLmDkyJEa20eOHIm4uLgWqpVpKygoAADY2toCAC5cuICysjKNGMrlcvj6+nIxjI+Ph6+vL+RyOVcmJCQESqWSa06Pj4/X+hxCQkLw559/oqysrN18Vq+//jpCQ0MxfPhwje0U58Zx4MABBAQEYPz48XBwcICfnx+2bt3K7b979y6ys7M17l8sFiM4OFgjzlKpFM899xxXJjAwEFKpVKNMcHAwxGIxVyYkJASZmZm4d+8eV0bXZ9EW4hwUFIQTJ07gxo0bAIBLly7h7NmzGD16NACKc1MxpbgaUhdDtMvkJDc3FxUVFZDJZBrbZTIZsrOzW6hWposxhoULFyIoKAi+vr4AgOzsbJiZmaFjx44aZavHMDs7WyvGHTt2hJmZWa1lZDIZysvLkZub2y4+qz179uDixYuIiorS2kdxbhx37tzB5s2b4eHhgaNHjyIsLAzz58/HN998AwDcPdZ2/9nZ2XBwcNA6t4ODQ51xrn4NfWXaQpwXL16MyZMnw9vbGyKRCH5+foiMjMTkyZMBUJybiinF1ZC6GMI05qltITweT+M9Y0xrGwHCw8Nx+fJlnD17ts6yNWOoK551lWF/d77i8XgaP9d2jtYqPT0dEREROHbsGCQSicHHUZzrR61WIyAgAKtXrwYA+Pn54erVq9i8eTOmT5/Olavr/hsaZ0Ov01rt3bsXO3fuxK5du9CzZ08kJSUhMjIScrkcM2bM4MpRnJuGKcW1obFvly0n9vb2EAgEWllcTk6OVrbX3r3xxhs4cOAATp48CRcXF267o6MjVCoV8vLyNMpXj6Gjo6NWjPPy8lBWVlZrmZycHAiFQtjZ2bX5z+rChQvIycmBv78/hEIhhEIhTp06hfXr10MoFEImk1GcG4GTkxN69Oihsc3HxwdpaWkAKuMDoNb7d3R0xIMHD7TO/fDhwzrjDKDOMm0hzosWLcKSJUswadIk9OrVC9OmTcOCBQu4VkGKc9MwpbgaUhdDtMvkxMzMDP7+/oiNjdXYHhsbiwEDBrRQrUwLYwzh4eHYt28ffv31V7i5uWns9/f3h0gk0ohhVlYWrly5wsWwf//+uHLlCrKysrgyx44dg1gshr+/P1em5udw7NgxBAQEQCQStfnPatiwYUhOTkZSUhL3CggIwJQpU7ifKc4NN3DgQK2h8Ddu3ICrqysAwM3NDY6Ojhr3r1KpcOrUKY04FxQUICEhgStz/vx5FBQUaJQ5ffq0xjDuY8eOQS6Xo2vXrlwZXZ9FW4hzSUkJ+HzNrxWBQMANJaY4Nw1TiqshdTGIwV1n25iqYZNfffUVS0lJYZGRkczS0pLdu3evpatmEv71r38xqVTKfvvtN5aVlcW9SkpKuDJhYWHMxcWFHT9+nF28eJENHTpU5xDXYcOGsYsXL7Ljx48zFxcXnUNcFyxYwFJSUthXX32ld4hre/msqo/WYYzi3BgSEhKYUChkq1atYjdv3mTffvsts7CwYDt37uTKREdHM6lUyvbt28eSk5PZ5MmTdQ7F7N27N4uPj2fx8fGsV69eGkMx8/PzmUwmY5MnT2bJycls3759zMbGRmMo5u+//84EAgGLjo5mqampLDo6ulUPca1uxowZzNnZmRtKvG/fPmZvb8/efvttrgzF2ThFRUUsMTGRJSYmMgBs7dq1LDExkRvqb0pxNaQudWm3yQljjG3cuJG5uroyMzMz9uyzz3LDZEnlUDVdr23btnFlSktLWXh4OLO1tWXm5uZszJgxLC0tTeM89+/fZ6Ghoczc3JzZ2tqy8PBwjWFqjDH222+/MT8/P2ZmZsa6du3KNm/erFWf9vRZ1UxOKM6N4+DBg8zX15eJxWLm7e3NtmzZorFfrVaz999/nzk6OjKxWMwGDx7MkpOTNco8evSITZkyhVlbWzNra2s2ZcoUlpeXp1Hm8uXLbNCgQUwsFjNHR0e2fPlybhhmle+//555eXkxkUjEvL29WUxMTNPcdDMrLCxkERERrEuXLkwikTB3d3e2bNkyplQquTIUZ+OcPHlS5+/kGTNmMMZMK66G1KUuPMZqTP1GCCGEENKC2mWfE0IIIYSYLkpOCCGEEGJSKDkhhBBCiEmh5IQQQgghJoWSE0IIIYSYFEpOCCGEEGJSKDkhhBBCiEmh5IQQQgghJoWSE0IIIYSYFEpOCCGEEGJSKDkhhBBCiEmh5IQQQgghJuX/A9dxe1yMX82GAAAAAElFTkSuQmCC",
      "text/plain": [
       "Figure(PyObject <Figure size 640x480 with 1 Axes>)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "6-element Array{PyCall.PyObject,1}:\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd00634150>\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd00634390>\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd00634550>\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd00634710>\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd006348d0>\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd00634a90>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plot(out[2])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "modelsol (generic function with 1 method)"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "function sicr!(du,u,p,t)\n",
    "\t# variables are concentrations X/N\n",
    "\tC,D,R,I,Z,U = u\n",
    "\tbeta,sigmac,sigmar,gammad,gammar = p\n",
    "\tdu[1] = dC = sigmac*I - gammar*C - gammad*C\n",
    "\tdu[2] = dD = gammad*C\n",
    "\tdu[3] = dR = gammar*C\n",
    "\tdu[4] = dI = beta*I*(1-Z) - sigmac*I - sigmar*I\n",
    "\tdu[5] = dZ = beta*I*(1-Z)\n",
    "\tdu[6] = dU = sigmar*I\n",
    "end\n",
    "\n",
    "function modelsol(pin,total,N)\n",
    "\tu0 = [0.,0.,0.,1/N,0.,0.]\n",
    "\tp = pin[1:end-1]\n",
    "\tif pin[end] < total - 1\n",
    "\t\tstart = ceil(Int,pin[end])\n",
    "\t\ttspan = (pin[end],total-1)\n",
    "\t\tprob = ODEProblem(sicr!,u0,tspan,p)\n",
    "\t\tsol = solve(prob,saveat=collect(start:total-1))\n",
    "\t\t# align solutions to time 0 of data and reorder rows and cols\n",
    "\t\tprediction = zeros(total,3)\n",
    "\t\tfor set in 1:3\n",
    "\t\t\tprediction[max(0,start)+1:end,set] .= sol[set,max(0,-start)+1:end]\n",
    "\t\tend\n",
    "\t\treturn prediction*N, sol\n",
    "\telse\n",
    "\t\treturn zeros(total,3)\n",
    "\tend\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "pin = ones(6)\n",
    "total = 60\n",
    "N = 6.67e7 #same as population\n",
    "\n",
    "pred, sol = modelsol(out[end],total,N);\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3df3RU5YHH/89lQoYfkll+SH5NKmhd/BG0bXAVdCStGrSiuDGrFZeV3R6OFIXEgHWx33Nke74lfoFiohRdf2ytrZq2yaBuVSR7inEooJBCBaTW1lCSMDGCmABCAsPz/SPmmiEzYfKDTG54v86ZY+Y+z3PzzAPO/XDvc59rGWOMAAAAHGRQvDsAAADQVQQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAnx7sCZcvLkSe3bt08jRoyQZVnx7g4AAIiBMUaHDh1SWlqaBg2Kfp5lwAaYffv2KSMjI97dAAAA3VBTUyOv1xu1fMAGmBEjRkhqHYCkpKQ49wYAAMSiqalJGRkZ9nE8mgEbYNouGyUlJRFgAABwmNNN/2ASLwAAcBwCDAAAcBwCDAAAcBwCDAAAcBwCDAAAcBwCDAAAcBwCDAAAcBwCDAAAcJwBu5AdAADofaFQSIFAQMFgUKmpqfL5fHK5XH3eDwIMAACIid/vV35+vmpra+1tXq9XJSUlys3N7dO+cAkJAACclt/vV15eXlh4kaS6ujrl5eXJ7/f3aX8IMAAAoFOhUEj5+fkyxnQoa9tWUFCgUCjUZ30iwAAAgE4FAoEOZ17aM8aopqZGgUCgz/pEgAEAAJ0KBoO9Wq83EGAAAECnUlNTe7VebyDAAACATvl8Pnm9XlmWFbHcsixlZGTI5/P1WZ8IMAAAoFMul0slJSWS1CHEtL0vLi7u0/VgCDAAAOC0cnNzVVZWpvT09LDtXq9XZWVlfb4OjGUi3RM1ADQ1Ncnj8aixsVFJSUnx7g4AAAPCmV6JN9bjNyvxAgCAmLlcLmVnZ8e7G1xCAgAAzkOAAQAAjkOAAQAAjkOAAQAAjkOAAQAAjkOAAQAAjkOAAQAAjkOAAQAAjkOAAQAAjkOAAQAAjtOlALNkyRJZlhX2SklJscuNMVqyZInS0tI0dOhQZWdna9euXWH7OHjwoGbNmiWPxyOPx6NZs2bp888/D6uzY8cOTZ06VUOHDlV6erp+/OMfa4A+sgkAAHRDl8/AXHrppQoGg/Zrx44ddtmyZcu0cuVKrVq1Slu2bFFKSopuuOEGHTp0yK4zc+ZMbd++XWvXrtXatWu1fft2zZo1yy5vamrSDTfcoLS0NG3ZskVPPPGEVqxYoZUrV/bwowIAgIGiyw9zTEhICDvr0sYYo+LiYv3oRz+yH6n9i1/8QsnJyXrppZd07733avfu3Vq7dq02b96sK6+8UpL0zDPPaPLkyfrwww81YcIEvfjiizp27Jief/55ud1uZWZm6i9/+YtWrlypwsJCWZbVw48MAACcrstnYD766COlpaVp/Pjx+t73vqePP/5YklRdXa36+nrl5OTYdd1ut6ZOnaqNGzdKkjZt2iSPx2OHF0m66qqr5PF4wupMnTpVbrfbrjNt2jTt27dPe/bsidqv5uZmNTU1hb0AAMDA1KUAc+WVV+qFF17QW2+9pWeeeUb19fWaMmWKDhw4oPr6eklScnJyWJvk5GS7rL6+XmPHju2w37Fjx4bVibSPtrJoioqK7Hk1Ho9HGRkZXfloAADAQboUYG666Sbdfvvtmjhxoq6//nq9/vrrklovFbU59RKPMSZsW6RLQKer0zaBt7PLR4sXL1ZjY6P9qqmp6cInAwAATtKj26iHDx+uiRMn6qOPPrLnxZx6lqShocE+g5KSkqJPPvmkw34+/fTTsDqR9iF1PLvTntvtVlJSUtgLAAAMTD0KMM3Nzdq9e7dSU1M1fvx4paSkqKKiwi5vaWlRZWWlpkyZIkmaPHmyGhsb9d5779l13n33XTU2NobVeeedd9TS0mLXWbdundLS0jRu3LiedBcAAAwQXQowixYtUmVlpaqrq/Xuu+8qLy9PTU1Nuueee2RZlgoKCrR06VKtWbNGO3fu1OzZszVs2DDNnDlTknTxxRfrxhtv1Jw5c7R582Zt3rxZc+bM0fTp0zVhwgRJrbdZu91uzZ49Wzt37tSaNWu0dOlS7kACAAC2Lt1GXVtbq7vuukv79+/Xueeeq6uuukqbN2/WeeedJ0n64Q9/qKNHj2revHk6ePCgrrzySq1bt04jRoyw9/Hiiy9qwYIF9t1Kt956q1atWmWXezweVVRU6L777tOkSZM0cuRIFRYWqrCwsDc+LwAAGAAsM0CXuG1qapLH41FjYyPzYQAAcIhYj988CwkAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADhOl9aBAQAAA18oFFIgEFAwGFRqaqp8Pp9cLle8uxWGAAMAAGx+v1/5+fmqra21t3m9XpWUlCg3NzeOPQvHJSQAACCpNbzk5eWFhRdJqqurU15envx+f5x61hEBBgAAKBQKKT8/X5EW6G/bVlBQoFAo1Nddi4gAAwAAFAgEOpx5ac8Yo5qaGgUCgT7sVXQEGAAAoGAw2Kv1zjQCDAAAUGpqaq/WO9MIMAAAQD6fT16vV5ZlRSy3LEsZGRny+Xx93LPICDAAAEAul0slJSWS1CHEtL0vLi7uN+vBEGAAAIAkKTc3V2VlZUpPTw/b7vV6VVZW1q/WgbFMpPulBoCmpiZ5PB41NjYqKSkp3t0BAMAx4rkSb6zHb1biBQAAYVwul7Kzs+PdjU5xCQkAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADhOjwJMUVGRLMtSQUGBva25uVnz58/XmDFjNHz4cN16662qra0Na7d3717dcsstGj58uMaMGaMFCxaopaUlrE5lZaWysrI0ZMgQnX/++Xrqqad60lUAADCAdDvAbNmyRU8//bQuu+yysO0FBQVas2aNSktLtWHDBh0+fFjTp09XKBSSJIVCId188806cuSINmzYoNLSUpWXl2vhwoX2Pqqrq/Xd735XPp9P27Zt08MPP6wFCxaovLy8u90FAAADiemGQ4cOmQsvvNBUVFSYqVOnmvz8fGOMMZ9//rkZPHiwKS0ttevW1dWZQYMGmbVr1xpjjHnjjTfMoEGDTF1dnV3n5ZdfNm632zQ2NhpjjPnhD39oLrroorDfee+995qrrroq5j42NjYaSfY+AQBA/xfr8btbZ2Duu+8+3Xzzzbr++uvDtldVVen48ePKycmxt6WlpSkzM1MbN26UJG3atEmZmZlKS0uz60ybNk3Nzc2qqqqy67TfR1udrVu36vjx4xH71NzcrKamprAXAAAYmLocYEpLS/XHP/5RRUVFHcrq6+uVmJiokSNHhm1PTk5WfX29XSc5OTmsfOTIkUpMTOy0TnJysk6cOKH9+/dH7FdRUZE8Ho/9ysjI6OpHAwAADtGlAFNTU6P8/Hz96le/0pAhQ2JuZ4yRZVn2+/Y/x1rHGBO1rSQtXrxYjY2N9qumpibm/gEAAGfpUoCpqqpSQ0ODsrKylJCQoISEBFVWVurxxx9XQkKCkpOT1dLSooMHD4a1a2hosM+opKSk2Gda2hw8eFDHjx/vtE5DQ4MSEhI0evToiH1zu91KSkoKewEAgIGpSwHmuuuu044dO7R9+3b7NWnSJN199932z4MHD1ZFRYXdJhgMaufOnZoyZYokafLkydq5c6eCwaBdZ926dXK73crKyrLrtN9HW522/QMAgLNbQlcqjxgxQpmZmWHbhg8frtGjR9vbv//972vhwoUaPXq0Ro0apUWLFmnixIn2hN+cnBxdcsklmjVrlpYvX67PPvtMixYt0pw5c+yzJnPnztWqVatUWFioOXPmaNOmTXruuef08ssv98ZnBgAADtelABOLxx57TAkJCbrjjjt09OhRXXfddXr++eflcrkkSS6XS6+//rrmzZunq6++WkOHDtXMmTO1YsUKex/jx4/XG2+8oQceeEA/+9nPlJaWpscff1y33357b3cXAAA4kGXaZscOME1NTfJ4PGpsbGQ+DAAADhHr8ZtnIQEAAMchwAAAAMchwAAAAMchwAAAAMfp9buQAABA/xYKhRQIBBQMBpWamiqfz2ffLewUBBgAAM4ifr9f+fn5qq2ttbd5vV6VlJQoNzc3jj3rGi4hAQBwlvD7/crLywsLL5JUV1envLw8+f3+OPWs6wgwAACcBUKhkPLz8xVp+be2bQUFBQqFQn3dtW4hwAAAcBYIBAIdzry0Z4xRTU2NAoFAH/aq+wgwAACcBdo/RLk36sUbAQYAgLNAampqr9aLNwIMAABnAZ/PJ6/XK8uyIpZblqWMjAz5fL4+7ln3EGAAADgLuFwulZSUSFKHENP2vri42DHrwRBgAAA4S+Tm5qqsrEzp6elh271er8rKyhy1DoxlIt1PNQDE+jhuAADONv15Jd5Yj9+sxAsAwFnG5XIpOzs73t3oES4hAQAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAxyHAAAAAx+lSgHnyySd12WWXKSkpSUlJSZo8ebLefPNNu7y5uVnz58/XmDFjNHz4cN16662qra0N28fevXt1yy23aPjw4RozZowWLFiglpaWsDqVlZXKysrSkCFDdP755+upp57qwUcEAAADTZcCjNfr1aOPPqqtW7dq69at+s53vqMZM2Zo165dkqSCggKtWbNGpaWl2rBhgw4fPqzp06crFApJkkKhkG6++WYdOXJEGzZsUGlpqcrLy7Vw4UL7d1RXV+u73/2ufD6ftm3bpocfflgLFixQeXl5L35sAADgZJYxxvRkB6NGjdLy5cuVl5enc889V7/85S915513SpL27dunjIwMvfHGG5o2bZrefPNNTZ8+XTU1NUpLS5MklZaWavbs2WpoaFBSUpIeeughvfbaa9q9e7f9O+bOnas//elP2rRpU8z9ampqksfjUWNjo5KSknryEQEAQB+J9fjd7TkwoVBIpaWlOnLkiCZPnqyqqiodP35cOTk5dp20tDRlZmZq48aNkqRNmzYpMzPTDi+SNG3aNDU3N6uqqsqu034fbXW2bt2q48ePR+1Pc3Ozmpqawl4AAGBg6nKA2bFjh8455xy53W7NnTtXa9as0SWXXKL6+nolJiZq5MiRYfWTk5NVX18vSaqvr1dycnJY+ciRI5WYmNhpneTkZJ04cUL79++P2q+ioiJ5PB77lZGR0dWPBgAAHKLLAWbChAnavn27Nm/erB/84Ae655579MEHH0Stb4yRZVn2+/Y/x1qn7SpXpLZtFi9erMbGRvtVU1MT82cCAADOktDVBomJifr6178uSZo0aZK2bNmikpIS3XnnnWppadHBgwfDzsI0NDRoypQpkqSUlBS9++67Yfs7ePCgjh8/bp91SUlJsc/GtN9HQkKCRo8eHbVfbrdbbre7qx8HAAA4UI/XgTHGqLm5WVlZWRo8eLAqKirssmAwqJ07d9oBZvLkydq5c6eCwaBdZ926dXK73crKyrLrtN9HW51JkyZp8ODBPe0uAAAYALp0Bubhhx/WTTfdpIyMDB06dEilpaV6++23tXbtWnk8Hn3/+9/XwoULNXr0aI0aNUqLFi3SxIkTdf3110uScnJydMkll2jWrFlavny5PvvsMy1atEhz5syxZxrPnTtXq1atUmFhoebMmaNNmzbpueee08svv9z7nx4AADhSlwLMJ598olmzZikYDMrj8eiyyy7T2rVrdcMNN0iSHnvsMSUkJOiOO+7Q0aNHdd111+n555+Xy+WSJLlcLr3++uuaN2+err76ag0dOlQzZ87UihUr7N8xfvx4vfHGG3rggQf0s5/9TGlpaXr88cd1++239+LHBgAATtbjdWD6K9aBAQDAec74OjAAAADxQoABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOQ4ABAACOkxDvDgAAgK4LhUIKBAIKBoNKTU2Vz+eTy+WKudzpCDAAADiM3+9Xfn6+amtr7W1er1clJSXKzc09bflAYBljTLw7cSY0NTXJ4/GosbFRSUlJ8e4OAAC9wu/3Ky8vT6cevi3LkiQtWrRIK1asiFpeVlbWr0NMrMdvAgwAAA4RCoU0bty4sDMrp3K5XAqFQhHLLMuS1+tVdXV1v72cFOvxm0m8AAA4RCAQ6DS8SIoaXiTJGKOamhoFAoHe7lqfI8AAAOAQwWCwX+0nnggwAAA4RGpqar/aTzwRYAAAcAifzyev12tPyI3E5XJFLbcsSxkZGfL5fGeqi32GAAMAgEO4XC6VlJRIUoeQYlmWLMtSYWFh1HJJKi4u7rcTeLuCAAMAgIPk5uaqrKxM6enpYdu9Xq/Kysq0bNmyTsv78y3UXcFt1AAAONBAXYmXdWAIMAAAOA7rwAAAgAGLAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAAByHAAMAABynSwGmqKhIV1xxhUaMGKGxY8fqtttu04cffhhWp7m5WfPnz9eYMWM0fPhw3XrrraqtrQ2rs3fvXt1yyy0aPny4xowZowULFqilpSWsTmVlpbKysjRkyBCdf/75euqpp7r5EQEAwEDTpQBTWVmp++67T5s3b1ZFRYVOnDihnJwcHTlyxK5TUFCgNWvWqLS0VBs2bNDhw4c1ffp0hUIhSVIoFNLNN9+sI0eOaMOGDSotLVV5ebkWLlxo76O6ulrf/e535fP5tG3bNj388MNasGCBysvLe+ljAwAAJ7OMMaa7jT/99FONHTtWlZWVuvbaa9XY2Khzzz1Xv/zlL3XnnXdKkvbt26eMjAy98cYbmjZtmt58801Nnz5dNTU1SktLkySVlpZq9uzZamhoUFJSkh566CG99tpr2r17t/275s6dqz/96U/atGlTTH1ramqSx+NRY2OjkpKSuvsRAQBAH4r1+N2jOTCNjY2SpFGjRkmSqqqqdPz4ceXk5Nh10tLSlJmZqY0bN0qSNm3apMzMTDu8SNK0adPU3Nysqqoqu077fbTV2bp1q44fPx6xL83NzWpqagp7AQCAganbAcYYo8LCQl1zzTXKzMyUJNXX1ysxMVEjR44Mq5ucnKz6+nq7TnJyclj5yJEjlZiY2Gmd5ORknThxQvv374/Yn6KiInk8HvuVkZHR3Y8GAAD6uW4HmPvvv1/vv/++Xn755dPWNcbIsiz7ffufY63TdqUrUltJWrx4sRobG+1XTU1NTJ8DAAA4T7cCzPz58/Xaa69p/fr18nq99vaUlBS1tLTo4MGDYfUbGhrsMyopKSn2mZY2Bw8e1PHjxzut09DQoISEBI0ePTpin9xut5KSksJeAABgYOpSgDHG6P7775ff79fvf/97jR8/Pqw8KytLgwcPVkVFhb0tGAxq586dmjJliiRp8uTJ2rlzp4LBoF1n3bp1crvdysrKsuu030dbnUmTJmnw4MFd+4QAAGDA6dJdSPPmzdNLL72kV199VRMmTLC3ezweDR06VJL0gx/8QL/73e/0/PPPa9SoUVq0aJEOHDigqqoquVwuhUIhfeMb31BycrKWL1+uzz77TLNnz9Ztt92mJ554QlLrbdSZmZm69957NWfOHG3atElz587Vyy+/rNtvvz2mvnIXEgAAzhPr8btLASba/JOf//znmj17tiTp2LFjevDBB/XSSy/p6NGjuu6667R69eqwSbV79+7VvHnz9Pvf/15Dhw7VzJkztWLFCrndbrtOZWWlHnjgAe3atUtpaWl66KGHNHfu3Fi7SoABADheKBRSIBBQMBhUamqqfD6fXC5XvLt1Rp2RAOMkBBgAgJP5/X7l5+eHrWbv9XpVUlKi3NzcOPbszOqTdWAAAEDv8/v9ysvL6/Aonrq6OuXl5cnv98epZ/0HAQYAgH4kFAopPz9fkS6QtG0rKCiwH9FztiLAAADQjwQCgQ5nXtozxqimpkaBQKAPe9X/EGAAAOhH2i8z0hv1BioCDAAA/Uhqamqv1huoCDAAAPQjPp9PXq836tIllmUpIyNDPp+vj3vWvxBgAADoR1wul0pKSiR1XH+t7X1xcfGAXw/mdAgwAAD0M7m5uSorK1N6enrYdq/Xq7KysgG9DkysWMgOAIB+ipV4ox+/E/qwTwAAoAtcLpeys7Pj3Y1+iUtIAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcQgwAADAcXgWEgAAcXI2PqyxtxBgAACIA7/fr/z8fNXW1trbvF6vSkpKlJubG8eeOQOXkAAA6GN+v195eXlh4UWS6urqlJeXJ7/fH6eeOQcBBgCAPhQKhZSfny9jTIeytm0FBQUKhUJ93TVHIcAAANCHAoFAhzMv7RljVFNTo0Ag0Ie9ch4CDAAAfSgYDPZqvbMVAQYAgD6Umpraq/XOVgQYAAD6kM/nk9frlWVZEcsty1JGRoZ8Pl8f98xZCDAAAPQhl8ulkpISSeoQYtreFxcXsx7MaRBgAADoY7m5uSorK1N6enrYdq/Xq7KyMtaBiYFlIt3HNQA0NTXJ4/GosbFRSUlJ8e4OAAAdsBJvR7Eev1mJFwCAOHG5XMrOzo53NxyJS0gAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDAAAMBxeBYSAAA9cLoHMvLAxjODAAMAQDf5/X7l5+ertrbW3ub1elVSUqLc3NzTlqP7LGOMiXcnzoRYH8cNAEB3+P1+5eXl6dTDqGVZkqRFixZpxYoVUcvLysoIMRHEevwmwAAA0EWhUEjjxo0LO7NyKpfLpVAoFLHMsix5vV5VV1dzOekUsR6/mcQLAEAXBQKBTsOLpKjhRZKMMaqpqVEgEOjtrp01CDAAAHRRMBjsV/s5GxFgAADootTU1H61n7NRlwPMO++8o1tuuUVpaWmyLEuvvPJKWLkxRkuWLFFaWpqGDh2q7Oxs7dq1K6zOwYMHNWvWLHk8Hnk8Hs2aNUuff/55WJ0dO3Zo6tSpGjp0qNLT0/XjH/+4w0QoAADiwefzyev12hNyI3G5XFHLLctSRkaGfD7fmerigNflAHPkyBFdfvnlWrVqVcTyZcuWaeXKlVq1apW2bNmilJQU3XDDDTp06JBdZ+bMmdq+fbvWrl2rtWvXavv27Zo1a5Zd3tTUpBtuuEFpaWnasmWLnnjiCa1YsUIrV67sxkcEAKB3uVwulZSUSFKHkGJZlizLUmFhYdRySSouLmYCb0+YHpBk1qxZY78/efKkSUlJMY8++qi97dixY8bj8ZinnnrKGGPMBx98YCSZzZs323U2bdpkJJk///nPxhhjVq9ebTwejzl27Jhdp6ioyKSlpZmTJ0/G1LfGxkYjyTQ2NvbkIwIAEFV5ebnxer1Gkv3KyMgw5eXlMZWjo1iP3706B6a6ulr19fXKycmxt7ndbk2dOlUbN26UJG3atEkej0dXXnmlXeeqq66Sx+MJqzN16lS53W67zrRp07Rv3z7t2bOnN7sMAEC35ebmas+ePVq/fr1eeuklrV+/XtXV1fb6LqcrR/f16kq89fX1kqTk5OSw7cnJyfr73/9u1xk7dmyHtmPHjrXb19fXa9y4cR320VY2fvz4Du2bm5vV3Nxsv29qaur+BwEAIEYul0vZ2dndLkf3nJG7kE693meMCdsWaVLT6eqYLyfwRpsQVVRUZE8K9ng8ysjI6Hb/AQBA/9arASYlJUXSV2di2jQ0NNhnUFJSUvTJJ590aPvpp5+G1Ym0D6nj2Z02ixcvVmNjo/2qqanp2YcBAAD9Vq8GmPHjxyslJUUVFRX2tpaWFlVWVmrKlCmSpMmTJ6uxsVHvvfeeXefdd99VY2NjWJ133nlHLS0tdp1169YpLS2tw6WlNm63W0lJSWEvAAAwMHU5wBw+fFjbt2/X9u3bJbVO3N2+fbv27t0ry7JUUFCgpUuXas2aNdq5c6dmz56tYcOGaebMmZKkiy++WDfeeKPmzJmjzZs3a/PmzZozZ46mT5+uCRMmSGq9zdrtdmv27NnauXOn1qxZo6VLl6qwsLDTe+4BAMBZoqu3N61fvz7sdrC21z333GOMab2V+pFHHjEpKSnG7Xaba6+91uzYsSNsHwcOHDB33323GTFihBkxYoS5++67zcGDB8PqvP/++8bn8xm3221SUlLMkiVLYr6F2hhuowYAwIliPX7zNGoAANBv8DRqAAAwYBFgAACA4xBgAACA4/TqSrwAAAxEoVBIgUBAwWBQqamp8vl8PIgxzggwAAB0wu/3Kz8/X7W1tfY2r9erkpISnmkUR1xCAgAgCr/fr7y8vLDwIkl1dXXKy8uT3++PU89AgAEAIIJQKNhmV6YAABZySURBVKT8/HxFWm2kbVtBQYFCoVBfdw0iwAAAEFEgEOhw5qU9Y4xqamoUCAT6sFdoQ4ABACCCYDDYq/XQuwgwAABEkJqa2qv10LsIMAAARODz+eT1eqM+RNiyLGVkZMjn8/VxzyARYAAAiMjlcqmkpESSOoSYtvfFxcWsBxMnBBgAAKLIzc1VWVmZ0tPTw7Z7vV6VlZWxDkwc8TRqAABOg5V4+06sx29W4gUA4DRcLpeys7Pj3Q20wyUkAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOAQYAADgOKzECwA4K3T2OAAeFeA8BBgAwIDn9/uVn5+v2tpae5vX67WfNh2tjIc19l88zBEAMKD5/X7l5eXp1MOdZVkdtrUvk8QTp+Mg1uM3c2AAAANWKBRSfn5+xKDS2b/f28oKCgoUCoXOWP/QfQQYAMCAFQgEwi4NdYUxRjU1NQoEAr3cK/QG5sAAAAaESBNxg8Fgj/fbG/tA7yPAAAAcL9ok3Tlz5vR436mpqT3eB3ofk3gBAI52ukm6o0eP1meffdbpnJdILMuS1+tVdXU1t1T3oViP35yBAQA4xqmXiaZMmdLpJN22u4mkjncdtX8fqUySiouLCS/9FJN4AQCO4Pf7NW7cOH3729/WzJkz9e1vf1vp6emdTtI1xujAgQNasmSJ0tPTw8q8Xq/Ky8tVXl4esYxbqCMLnQzp7T1v6+UdL+vtPW8rdDI+d2lxBgYA0O9Fu0y0f//+mNpfeOGF2rNnT9TVdmfMmMFKvDHw7/Yrf22+apvazTVK8qrkxhLlXty3YY85MACAfiPSnUSSNG7cuG7fDi1J69evV3Z2di/18uzk3+1X3m/yZHTKXCN9uejfHWW9EmKYAwMAcJTO7iTqbnhpm4jbFoQQm9DJkAJ7AwoeCip1RKqmeKcof21+h/AiSUZGliwVrC3QjAkz5BrUN2euCDAAgLiLdomorq5OjzzySLf2yUTc7ol0mWjMsDHa/0X0y3VGRjVNNQrsDSh7XHYf9JJJvACAOOvucv+nOvfcc8PeMxG369ouE7UPL5I6DS/tBQ/13aJ/nIEBAPSpU+e5hEKhHs1vabtM9Ne//lUbN25kIm43hU6Gol4milXqiL5b9I8AAwDoM5HmuYwaNSrm9p2t15KYmMhE3RidOsfF9zWfAnsDHc68xMqSJW+SV76v9d1cIwIMAKBXRbqTyOVyRZ3n8tlnn8W03//6r//SM88802GSb3FxMZeJuiDardB5l+R1a39tdyEV31jcZxN4JW6jBgD0omh3Eq1cuVKFhYXdulTUfkl/SazXEqNIZ1le/fDVqLdCx3rp6Nxh5+rTLz6132ckZaj4xuJeWweG26gBAH2qszuJ7rjjjm7tM9KdRFwmOr1IZ1nSR6Tr2IljUW+FliSX5dJJczJinbbLRH+d/1dtrN0YFoz68sxLGwIMAKDLuvpMoliNGjUq7JISl4i6LtqCc3WH6k7bNmRaHwtw6hmZ9peJEhMS++xW6c4QYAAAXRLpMtGYMWNiXta/M7/5zW/kcrm4RNRNvXEnUcGVBSrbXdZhjkxvXibqDQQYAEAHXZ2I29Pw0jbPJTs7m8ASo96+k6jNjItmaEXOig77jsdlos4QYAAAYU43Eben9350div02RZeIoWQ9kEhWrl/t1/5b+ar9lC7P6MRXuVd2r07iaTwW6FdRsreIykoKVVSRvtOh6RAQAoGpdRUyeeT4vDnxl1IAHCWinSW5dVXX414huXU0NFVbWdYVq5cqQceeCAsHGVkZJyV81xO92TnaCHlrol3acXG5TJG+nJqiiTJMpJp975TEdrKksruKFfubkn5+VL7O8a8XqmkpPXnaGW99OcX6/GbAAMADhXtMk8s5ZHOsqSnp+vYsWM6cOBAr/az7QxL27L+p+v32eB0T3ZeNGVRxJASVj1SWDGSS5ZOykQMM5aRRn0hDT0h1Xq+2p7RKBWvlXKnPyitWCGdGg0sq+O29mWSVFbWKyGGAEOAQT/SkwNNvNqe9jRxZ+W0DWsbamnRjtWr9cXf/qZhF1ygifPmyZWYeNqyzsr9fr8eWLBA4+vqlKrWM/3V6el67PHHlZub22m5JOXl5ckyRj7JLg9IOvnl7x0kRS07XXnymDG6aP9+u2yP16uVJSVfnWGJ15/D8RYFXl+t4Cd/U2ryBfLdPE+uwYmnLevttlNuulcXrP7HTuequDRIIXMyckiJhWltak49yyKp7DfSjD9LgfOk4DlS6mHJ93fJZdQ6XqFQ13+fZbWeiamu7vHlJALMGQgwLceO6oVVD6ou+JHSUy/Uv92/XIlDhsZU7sS2/bVfZ7LtkUNN+v/+n3/VSR3WuIxLeuX3+v1+FeTP16iEfRp2jvTFYemzE2kqLnnCPtBEK5cUl7a5kk4ULNAGV539BXdNKF0JxY+3/gvL749eTtuwtpt/+EOlP/ZT/c170i6/oHaQ6h5YKElRy65atixq2/U336JXX31VP7WkPe0OQuP+Li000gUPPqi/LV8esbzQSO+MHq1rDxzQyijlkqKWrZH0z1HKFxpp9OjRWj3UrQ0J+6KOZahggQLtxssXSpfry7GMWtbDtv6f/1D5H6xU7TlfHZy9h10quaRQkqKW5f77sl5vO+boIO0f2j4O9r6CjVLZpVHOsuw+g794/Xqph+v0DIgAs3r1ai1fvlzBYFCXXnqpiouL5fPF9pyF3g4wjy6+TY83v6pgu78MqY3SAvcM/WfRK52WS3Jc2/7arzPZ9qG5U/XLYe/06u/NO3CF3tuyRXtvVIeyr62Vrp3+oN753fKI5RlrW3+u6eO216+VZkgquDH8y8/75Zff7dMfVPnvlkctF23tssx/nKEdf3lVD0Qof2xt67+QI7V9bK005pwrtP/wlm61rf6zpfEXmaht/7ZbuuBidWvf3W1b8uXlCf/vlis/SrmkqGU9aasZM5TnfrX16supc0bUToSzFYuOX6EVg7f0attT55+cCeufbz2rEvEsy5n00kvSXXf1aBeODzC//vWvNWvWLK1evVpXX321/vu//1vPPvusPvjgA33ta187bfveDDCPLr5ND0f5yy9Jd3ycrN+c/0m3/oL3x7b9tV+na/svf0vWby/oXttp24fqrW8c7VbbWz706H8nNHar7bV/kN65Wo5qu/AP0k9pG1Pbwj9IK/tZ23j2a9EmSysmmz5taySNPiodGBpeZjMd27UvcxkpZJ2BtmeIZSRvk1Rd3AdhJRLOwEhXXnmlvvWtb+nJJ5+0t1188cW67bbbVFRUdNr2vRVgWo4d1biHhymYpL7/Cx7Htv21X7SlLW173ra/9utMte3LANErvvy8J60odxV9+XktRZ7j8tvfSrkfdG+oNMiSTpquD1kc5sAM6tFvOUNaWlpUVVWlnJycsO05OTnauHFjxDbNzc1qamoKe/WGF1Y92HqKPdqfpiWFBnVe7sS2/bVftKUtbXvetr/260y17ddOOYVgmdYuF2786n2k8kV/kNJPOcx5m6Tf/Eb66z+07vbUWTYn2/26aGXLrjIR29p3Gp36c/v3xcV9uh5Mvwww+/fvVygUUnJyctj25ORk1dfXR2xTVFQkj8djvzIyMiLW66q64Ee9sh8AAE517hfh771NrXcJLfs/6be/iR5SvtYo/eHZ1rkuL5W1/jfwnPS7f5RWXC39x4zWuS/t7RvRuv3fOymL1lZer1Re3vpKT+9Y1ku3UHdFv16J1zol5RljOmxrs3jxYhUWFtrvm5qaeiXEpKdeKB1Z1+P9AADiKNqlpBguXZ3uUk532rbNVflribTxa6dMtJUlZXh1+49X6rbCgvA7q0565Xrsy0XjQiF97ZTbxp/v5JZzb/vyzsoilIfdkj5jBivxRtPS0qJhw4bpt7/9rf75n//Z3p6fn6/t27ersrLytPvo7Tkw9UnR/wKfqb/g8WzbX/tF2/i3bVuxk7Zntu0g09quW3++ir6QWSzlPdl3f2trGWmUGaLPrGOtVa3wsvZNI80nabuTqLfbSq1nWnL/bClsgbhTF4XrJ8v29yVHz4FJTExUVlaWKioqwrZXVFRoypQpfduXIUPt22ijXYvMq06OWh7p5/7etr/2ayC2LTxxhaPaStKkE1e0/mOWtqdte+vYGd1uu/CaBztva0X587WkwqsXdbu8bT7JQGkrS3r6ey+q7LwHlX4k/MDvPeJS+XkPqjxKWdl5D2rZT947I23LzntQuf9vDJdjXK7Wu3ruuqv1vwM8vHRFvwwwklRYWKhnn31W//M//6Pdu3frgQce0N69ezV37tw+78t/Fr2ipc0zlHLKtcjUJmlp8wyVvlAftbyoeYaKHNa2v/bLqW3j8eXY07YPnvegBh0J/3oYdMSlB897UO/95L1Oy2n7Vdkr972i8jvLlZ4UfpBKT/Kq/M7yTsuW3bCs87Z3lCs9yRtW5vV4VXZHa9uybpaX39H9fffXtrkX5yr335dpz9IvtP7yx/RSyv1af/ljql76hXL/fVmnZZLOWFvl5kp79rTeevzSS63/ra7u87kkTtUvLyG1Wb16tZYtW6ZgMKjMzEw99thjuvbaa2Nqy0q8rMTbn9r25TLlvdW25XiLVr++Wn/75G+6IPkCzbt5nhLbte2snLbhbTt74nB3n0bc07b9tV9nsi2cwfHrwPQUz0ICAMB5HD0HBgAAoDMEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgEGAAA4DgJ8e7AmdL2jMqmpqY49wQAAMSq7bh9umdND9gAc+jQIUlSRkZGnHsCAAC66tChQ/J4PFHLLXO6iONQJ0+e1L59+zRixAhZltVr+21qalJGRoZqamo6fcw3WjFesWOsYsdYxY6xih1jFbszOVbGGB06dEhpaWkaNCj6TJcBewZm0KBB8nq9Z2z/SUlJ/AXvAsYrdoxV7Bir2DFWsWOsYnemxqqzMy9tmMQLAAAchwADAAAcx7VkyZIl8e6E07hcLmVnZyshYcBegetVjFfsGKvYMVaxY6xix1jFLt5jNWAn8QIAgIGLS0gAAMBxCDAAAMBxCDAAAMBxCDAAAMBxCDBdtHr1ao0fP15DhgxRVlaWAoFAvLsUd++8845uueUWpaWlybIsvfLKK2HlxhgtWbJEaWlpGjp0qLKzs7Vr16449Ta+ioqKdMUVV2jEiBEaO3asbrvtNn344YdhdZqbmzV//nyNGTNGw4cP16233qra2to49Th+nnzySV122WX2QlmTJ0/Wm2++aZczTtEVFRXJsiwVFBTY2xivVkuWLJFlWWGvlJQUu5zvq3B1dXX613/9V40ePVrDhg3TN77xDVVVVdnl8RwvAkwX/PrXv1ZBQYF+9KMfadu2bfL5fLrpppu0d+/eeHctro4cOaLLL79cq1atili+bNkyrVy5UqtWrdKWLVuUkpKiG264wX5e1dmksrJS9913nzZv3qyKigqdOHFCOTk5OnLkiF2noKBAa9asUWlpqTZs2KDDhw9r+vTpCoVCcex53/N6vXr00Ue1detWbd26Vd/5znc0Y8YM+8uRcYpsy5Ytevrpp3XZZZeFbWe8vnLppZcqGAzarx07dthlfF995eDBg7r66qs1ePBgvfnmm/rggw/005/+VP/wD/9g14nreBnE7J/+6Z/M3Llzw7ZddNFF5j//8z/j1KP+R5JZs2aN/f7kyZMmJSXFPProo/a2Y8eOGY/HY5566ql4dLFfaWhoMJJMZWWlMcaYzz//3AwePNiUlpbaderq6sygQYPM2rVr49XNfmPkyJHm2WefZZyiOHTokLnwwgtNRUWFmTp1qsnPzzfG8PeqvUceecRcfvnlEcv4vgr30EMPmWuuuSZqebzHizMwMWppaVFVVZVycnLCtufk5Gjjxo1x6lX/V11drfr6+rBxc7vdmjp1KuMmqbGxUZI0atQoSVJVVZWOHz8eNl5paWnKzMw8q8crFAqptLRUR44c0eTJkxmnKO677z7dfPPNuv7668O2M17hPvroI6WlpWn8+PH63ve+p48//lgS31eneu211zRp0iT9y7/8i8aOHatvfvObeuaZZ+zyeI8XASZG+/fvVygUUnJyctj25ORk1dfXx6lX/V/b2DBuHRljVFhYqGuuuUaZmZmSWscrMTFRI0eODKt7to7Xjh07dM4558jtdmvu3Llas2aNLrnkEsYpgtLSUv3xj39UUVFRhzLG6ytXXnmlXnjhBb311lt65plnVF9frylTpujAgQN8X53i448/1pNPPqkLL7xQb731lubOnasFCxbohRdekBT/73fWSu4iy7LC3htjOmxDR4xbR/fff7/ef/99bdiw4bR1z9bxmjBhgrZv367PP/9c5eXluueee1RZWRm1/tk6TjU1NcrPz9e6des0ZMiQmNudjeN100032T9PnDhRkydP1gUXXKBf/OIXuuqqqyTxfdXm5MmTmjRpkpYuXSpJ+uY3v6ldu3bpySef1L/927/Z9eI1XpyBidGYMWPkcrk6pMqGhoYO6RNfaZvdz7iFmz9/vl577TWtX79eXq/X3p6SkqKWlhYdPHgwrP7ZOl6JiYn6+te/rkmTJqmoqEiXX365SkpKGKdTVFVVqaGhQVlZWUpISFBCQoIqKyv1+OOPKyEhQcnJyYxXFMOHD9fEiRP10Ucf8X11itTUVF1yySVh2y6++GL7xpV4jxcBJkaJiYnKyspSRUVF2PaKigpNmTIlTr3q/8aPH6+UlJSwcWtpaVFlZeVZOW7GGN1///3y+/36/e9/r/Hjx4eVZ2VlafDgwWHjFQwGtXPnzrNyvE5ljFFzczPjdIrrrrtOO3bs0Pbt2+3XpEmTdPfdd9s/M16RNTc3a/fu3UpNTeX76hRXX311h2Ue/vKXv+i8886T1A++38/4NOEBpLS01AwePNg899xz5oMPPjAFBQVm+PDhZs+ePfHuWlwdOnTIbNu2zWzbts1IMitXrjTbtm0zf//7340xxjz66KPG4/EYv99vduzYYe666y6Tmppqmpqa4tzzvveDH/zAeDwe8/bbb5tgMGi/vvjiC7vO3LlzjdfrNf/3f/9n/vjHP5rvfOc75vLLLzcnTpyIY8/73uLFi80777xjqqurzfvvv28efvhhM2jQILNu3TpjDON0Ou3vQjKG8WqzcOFC8/bbb5uPP/7YbN682UyfPt2MGDHC/h7n++or7733nklISDA/+clPzEcffWRefPFFM2zYMPOrX/3KrhPP8SLAdNHPfvYzc95555nExETzrW99y7799Wy2fv16I6nD65577jHGtN5q98gjj5iUlBTjdrvNtddea3bs2BHfTsdJpHGSZH7+85/bdY4ePWruv/9+M2rUKDN06FAzffp0s3fv3vh1Ok7+4z/+w/5/7dxzzzXXXXedHV6MYZxO59QAw3i1uvPOO01qaqoZPHiwSUtLM7m5uWbXrl12Od9X4f73f//XZGZmGrfbbS666CLz9NNPh5XHc7wsY4w58+d5AAAAeg9zYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOMQYAAAgOP8/0ZvmFz5hp7KAAAAAElFTkSuQmCC",
      "text/plain": [
       "Figure(PyObject <Figure size 640x480 with 1 Axes>)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "1-element Array{PyCall.PyObject,1}:\n",
       " PyObject <matplotlib.lines.Line2D object at 0x7fbd52716150>"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plot(pred[:,1],color=\"black\")\n",
    "print(pred([:,1]))\n",
    "plot(pred[:,2],color=\"red\")\n",
    "plot(pred[:,3],color=\"green\")\n",
    "plot(df[:,:active],color=\"black\",marker=\"o\",linestyle=\"none\")\n",
    "plot(df[:,:cum_recover],color=\"red\",marker=\"o\",linestyle=\"none\")\n",
    "plot(df[:,:cum_deaths],color=\"green\",marker=\"o\",linestyle=\"none\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 1.4.0",
   "language": "julia",
   "name": "julia-1.4"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "1.4.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
